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Abstract. We have compiled a sample of 165 radio galax- 
ies from the literature to study the properties of the ex- 
tended emission line regions and their interaction with the 
radio source over a large range of redshift < z < 5.2. For 
each source, we have collected radio (size, lobe distance ra- 
tio and power) and spectroscopic parameters (luminosity, 
line width and equivalent width) for the four brightest UV 
lines. We also introduce a parameter Aj^ya measuring the 
asymmetry of the Lya line, assuming the intrinsic redshift 
of the line is the same as that for the He II A 1640 line, 
and show that this parameter is a good measure of the 
amount of absorption in the Lya line. 

Using these 18 parameters, we examine the statisti- 
cal significance of all 153 mutual correlations, and find 
the following significant correlations: (i) Lya asymmetry 
^LyQwith radio size D and redshift z, (ii) line luminosity 
with radio power, (iii) line luminosities of Lya, C IV, He II 
and C III] with each other, and (iv) equivalent widths of 
Lya, C IV, He II and C III] with each other. We interpret 
the correlation between redshift and Aj^ya as an increase in 
the amount of Hi around radio galaxies at 2; > 3. The al- 
most exclusive occurrence of H i absorption in small radio 
sources could indicate a denser surrounding medium or an 
un-pressurized, low density region, as suggested by Binette 
et al. (2000). Correlations (ii) to (iv) provide evidence for 
a common energy source for the radio power and total 
emission line luminosity, as found in flux density-limited 
samples of radio sources. 

The luminosity of the Lya line relative to the other 
emission lines and the continuum shows a strong increase 
at z > 3, coincident with the increase in the amount of as- 
sociated H I absorption. This indicates an increased abun- 
dance of hydrogen, both ionized and neutral, which may 
well be the reservoir of primordial hydrogen from which 
the galaxy is forming. This metallicity evolution is also 
seen in the nitrogen abundance, which shows a variation 
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of more than an order of magnitude, with the z > 3 radio 
galaxies occupying only the Z < 2Zq region. 

To examine the ionization mechanism of the extended 
emission line regions in HzRGs, we plot the UV emission 
line data in line-ratio diagnostic diagrams. The diagrams 
involving the high ionization C IV, He II and C III] lines 
seem to confirm previous results showing that AGN photo- 
ionization provides the best fit to the data. However, these 
models cannot fit the C II]/C III] ratio, which lies closer 
to the predictions for the highest velocity shock ioniza- 
tion models. We note that the C II] line is five times 
more sensitive to shock ionization than the high ioniza- 
tion UV lines, and show that a combination of shock and 
photo-ionization provides a better overall fit to the inte- 
grated spectra of HzRGs. A substantial contribution from 
shock ionization will show up first in shock sensitive lines 
like C II] or Mg II. We also confirm the findings of Best, 
Rottgering & Longair (2000b) that shock ionization occurs 
almost exclusively in small radio sources, and show that 
the angular size distribution can indeed explain the dif- 
ferences in three HzRG composite spectra. Because most 
HzRGs have radio sizes < 150 kpc, their integrated spectra 
might well contain a significant contribution from shock 
ionized emission. 
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1. Introduction 

During the last two decades, high redshift radio galaxies 
(HzRGs) have been used as probes of galaxy for- 
mation. Out to z ^ 1, they are uniquely identified 
with massive ellipticals (e.g., Lilly & Longair 1984; 



Best, Longair & Rottgering 1998). There are strong 



indications that this is also true at higher redshifts, 
mostly based on the remarkably strong correla- 
tion in the Hubble K — z diagram out to z = 5.19 



http: / /link. springer.de/link / service/journals / 00230 / index.htm 



(Lilly 1989; 


Eales et al. 1997 




van Breugel et al. 1998 


v^an Breugel et al. 199!;). By studying radio galaxies over 
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a large range of redshift, we can thus study the formation 
and evolution of massive galaxies. 

The spectra of HzRGs are mostly dominated by char- 
acteristic extended emission lines, indicative of a halo of 
ionized gas. The most prominent line in these extended 
emission line regions is Lya: the luminosity can reach ~ 
10''^ e rg s~^ and the spatial extent can be up to ~ 150 kpc 
(e.g., [van Ojik et al. 1996i|Adam et al. 1997| ). This allows 



a detailed study of the kinematics of the gas, both close 
to the galaxy, where interactions with the radio jets can 
be studied, as well as at large distances from the AGN 
where the gas is still undisturbed and should trace the 
primordial distribution of the gas in these massive galax- 
ies. Within the extent of the radio source, the high line 
velocity widths 1500 km s^^) and disturbed morpholo- 



gies (e.g., Villar-Martm et al. 1999b, Bicknell et al. 2000) 



of the emission line gas indicate a strong interaction with 
the radio jets. This interaction is also obvious in statistical 
comparisons of the radio and optical morphologies in z > 
0.6 radio galaxies: (i) the UV and optical emission is often 



aligned with the radio source (e.g., 


McCarthy et al. 1987, 


Chambers, Miley & van Breugel 1987), (ii) the radio and 



emission-line gas morphological asymmetries are strongly 



correlated (McCarthy, van Breugel & Kapahi 1991a), and 



(iii) the rest-f rame U— or i?— band morpholog y depends 
on radio size ( Best, Longair fc Rottgcring 1996| ). 

The bright Lya emission profiles frequently show 
narrow absorption caused by neutral Hi in the HzRG 



(van Ojik et al. 1997). This Hi absorption seems to occur 
almost exclusively in small radio sources, prompting van 
Ojik et al. to suggest that the H i absorption indicates a 
dense intergalactic region which confines the radio source. 
However, the detection of highly ionized C IV AA 1549 
absorption in 0943—242 at z = 2.93 suggests that the 
absorption is located further out in a low-metallicity gas 



shell (Binctte ct al. 2000), providing us with a new tool 
to probe the outer regions of forming massive galaxies. 

The radio power and emission line luminosities 
are found to be correlated (e 



[McCarthy 1993 

Willott ct al. 1999| ), indicating that the central AGN is 



the common energy source for both. The most likely 
mechanisms for the transformation of this AGN en- 
ergy into emission line luminosity are photo-ionization 
by an anisotropic UV source and shock excitation. Be- 
cause the ionizing spectra of these mechanisms are 
quite different, they will lead to differences in the 
emission line spectra, which can be used to deter- 
mine the dominating source of ionization. Such stud- 
ies using UV line ratio diagrams suggest that the 
main ionizing mechanism in HzRGs is nuclear photo- 



ionization (e.g.. 


Villar-Marti'n, Tadhuntcr & Clark 1997, 


Allen, Dopita & Tsvetanov 1998 


, hereafter VM97 and 



1 ra- 



ADT98). However, a recent study of z 
dio galaxies finds that shock ionization is impor- 
t ant in most of the smaller (< 150 kpc) sources 
([Best, Rottgcring fc Longair 2000bt hereafter BRLOO). 



Because the radio sources in HzRGs generally have sizes 
< 150 kpc, this is in apparent contradiction with the re- 
sults of VM97 and ADT98 (unless there is a drastic change 
in the ionization mechanism between z ^ 1 and z > 2). 

From the above, it is clear that the properties 
of the gas in HzRGs can provide an important di- 
agnostic to study various processes in forming mas- 
sive galaxies, such as star formation, chemical enrich- 
ment of the interstellar matter, and the influence of the 
AGN through shocks or photo-ionization. Detailed ob- 
servations of several representative individual galaxies 



(such as 4C 41.17, 


Chambcrs, Miley & van Breugel 1990 


pey et al. 1997, 


Bicknell et al. 2000|) can be used to deter- 



mine the relative importance of these mechanisms, but 
these should be complemented with a search for statistical 
relations between the emission line properties and other 
radio properties of a large sample of HzRGs. Such a study 
for z < 2.5 radio galaxies ( Baum fc McCarthy 200C| ) finds 



systematically larger line widths and velocity field ampli- 
tudes at z > 0.6 than at lower redshifts, but remains in- 
conclusive on the origin: gravitational or due to jet-cloud 
interactions. 

To allow such a statistical study over a continuous red- 
shift range z = 0toz = 5.2, we have compiled from the 
literature the emission line properties of a large sample of 
radio galaxies, together with relevant radio properties. In 
this paper, we first describe the compilation of our HzRG 
sample (§2), and then go on to discuss the determination 
of the different radio and spectroscopic parameters (§3). 
In §4, we perform a statistical analysis of correlations be- 
tween these parameters, taking the sometimes strong se- 
lection effects into account. In §5, we use diagnostic line- 
ratio diagrams to examine the ionization mechanisms in 
HzRG, and we discuss the implications of our results for 
the nature of HzRGs in §6. We present our conclusions in 
§7. 

Throughout, we shall assume Hq — 50 km s ^Mpc ^, 
(7o=0.5, and A = 0, unless otherwise stated, but our 
results do not depend on the adopted cosmology. We 
shall abbreviate the emission lines as follows: N V for 
N V A 1240, C IV for C IV AA 1549, He II for He II A 1640, 
C III] for C III] A 1909, C II] for C II] A 2326, Mg II for 
Mg II A 2800, [O II] for [O II] A 3727, and [O III] for 
[O III] A 5007. 

2. Sample selection 

Table ^ lists nine samples designed to find 
HzRGs. The surveys can be divided into two 
classes: (i) several large flux density-limited sur- 
veys such as the 3CR (Spinrad et al. 1985), BRL 
s ample ([Best, Rottgering fc Lehnert 1999D , MRC 
( [McCarthy et al. 1996|), and 6C, 7C and 8C sur- 
veys (e.g., Lacy et al. 1999); (ii) several "filtered" 
surveys, which have been designed to find the 
highest redshift objects. For the latter, the radio 
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spectral index is most often used (ultra steep spec- 
trum sources, e.g., De Breuck et al. 2000a), some- 
times in combination with an angular size upper 
limit (e.g., Blundell et al. 1998). Alternatively, the 
filter consists of a flux density interval centered 
around the peak in the source counts around 1 Jy 



(e.g., [AUington-Smith 19821 ). 

These two types of samples are complementary, in the 
sense that the complete surveys can provide information 
on the entire range of values of the parameters that were 
used as a high redshift filter (usually spectral index or ra- 
dio size), while the filtered samples can extend the radio 
power and redshift coverage of the un-filtered samples. 
For example, the addition of data from filtered surveys 
partially compensates for the Malmquist bias in the flux 
density-limited surveys (see Fig. |l|). We therefore compiled 
spectroscopic data on HzRGs from the samples listed in 
Table ||, augmented with four sources from other small 
surveys. We shall concentrate on radio galaxies at 2; > 2, 
because above this redshift, the brightest UV lines (Lya, 
C IV, He H, and C HI]) can be observed with ground 
based optical spectrographs. 

Of the 145 known z > 2 radio galaxies, 78 have 
published spectroscopic parameters including line fluxes, 
equivalent widths, and line widths. In order to fully inves- 
tigate redshift dependence, we consider also spectroscopic 
data for the z < 2 sources from the samples that provided 
sources at z > 2. Our final sample contains a similar num- 
ber of sources at z < 2 and at z > 2 (Fig. ||). From Figure 
|l|, we can see that we indeed include sources at z > 2 with 
radio powers that are more than an order of magnitude 
lower than in the un-filtered samples. In appendix A, we 
list all 167 sources. 



3. Determination of source parameters 

We now describe the radio and spectroscopic parameters 
from the literature papers and some derived quantities, as 
listed in Table A.l. 



3.1. Radio parameters 
3.1.1. Radio power 

In order to obtain a uniform measurement of the ra- 
dio flux densities, we determined the low frequency ra- 
dio flux density from the 325 MHz WENSS {5 > +2S°; 
Rengclink ct al. 1997) or the 365 MHz Texas radio survey 
(-35°< 5 < +28°; pouglas et al. 1996|), and the 1.4 GH z 
radio flux density from the NVSS ( [Condon et al. 199^ . 
This procedure flnds radio flux densities for 88% of the 
sources in our HzRG sample. Of the 19 missing sources, 
one (USS 0529-549) lies outside the area covered by the 
surveys and one (VLA J123642-I-621331) is too faint to 
be detected in any of the three surveys. The remaining 
sources have not been detected in the incomplete Texas 



37 



36 



X 

'[/) 35 



S> 34 
o 

33 
32 



Redshift z 

Fig. 1. Radio power at 325 MHz against redshift. Sources 
from un-filtered surveys are plotted with -I- signs, and 
sources from samples with spectral index and/or radio size 
filters are plotted as filled circles. Note the tight correla- 
tion at z < 1 due to the fiux density limited 3C survey, 
and the filtered surveys that fill up the lower power regions 
at 1 < z < 4. 



catalogue. Because no other deep large area radio surveys 
at frequencies below 325 MHz or above 1.4 GHz existQ 
we could not correct for the spectral curvature when per- 
forming a K-correction in the calculation of the rest-frame 
mono-chromatic radio power. Because the radio spectra of 
HzRGs are predominantly concave, neglecting this spec- 
tral curvature will generally lead to an over-estimation of 
the radio power of the highest redshift sources. Neglect- 
ing a spectral curvature of Aa — 0.4 between rest and 
observed frequency at z = 5 will overestimate the radio 
power up to a factor of two. From the spectral curvature in 
a large USS sample (Figure 9 of De Breuck et al. 2000a) , 
we expect only 30% of the sources in our radio galaxy sam- 
ple to have more spectral curvature, while all our sources 
are at lower redshifts. We therefore believe that our rest- 
frame radio powers are accurate to within a factor of two. 

3.1.2. Radio size 

To compile the data on the projected radio sizes, we gen- 
erally took the largest separation between the components 
of the radio source from the same reference as the spec- 
troscopy data (listed in appendix A). For the sources of 



^ The Cambridge 6C, 7C and 8C surveys only cover 
areas of the sky at 5 >30°, 20° and 60°, respectively 



(Hales, Baldwin fc Warner 1993 ^ Riley, Waldram fc Riley 1999, Rees 1990), 



while the Texas surveys covers 5 > —35" (Douglas et al. 1996 J 



4 Carlos De Breuck et al.: Statistical study of emission lines from HzRGs 



Table 1. Radio surveys used to construct our sample of radio galaxies. 



Survey 


Flux density limit 


Spectral index limit 


Angular size limit 




Reference 


3C 


5i78 > 10 Jy 


none 


none 


15 


Laing, Riley & Longair 1983| 


BRL 


5*408 > 5 Jy 


none 


none 


39 


Best, RottgcriuK & Lehnert 1999| 


MRC 


S408 > 0.95 Jy 


none 


none 


19 


McCarthy et al. 1996 


B2 1 Jy 


1 < 5408 < 2 Jy 


none 


none 


3 


AlliuKton-Smith 1982 


6C/7C/8C 


varies" 


none 


none 


7 


Lacy et al. 1998| 


6C* 


0.96 < 5i5i < 2.00 Jy 


atlT < -0.981 


e < 15" 


2 


Blundell et al. 1998 


WN/TN 


5i4oo > 10 mjy 


^1400 , -1 on 
"325 < -l.oO 


e < 1' 


34 


De Breuck et al. 2000a| 


uss 


varies*" 


a < -1.0* 


none 


30 


RottKeriuK et al. 1994| 


MG 


55000 > 50 mJy 


4800 , n 

«i4oo < -0.75 


e < 10" 


14 


Stern et al. 1999| 



* The Cambridge surveys listed in Table 1 of Lacy et al. (1999) consist of five sub-samples, each having different flux density 
limits. 

'' The USS sample of Rottgering et al. (1994) consists of nine sub-samples, each with different flux density and spectral index 
limits. See Table 4 of Rottgering et al. (1994) for details. 
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Fig. 2. Redshift distribution of the sample of z > 2 radio 
galaxies with spectroscopic data, and the lower redshift 
extension. 



Rottgering et al. (1997), we used the radio sizes from 
Rottgering et al. (1994), for the MRC/1 Jy sample of Mc- 
Carthy et al. (1996) the sizes are from Kapahi et al. (1998), 
and for the MG sample of Stern et al. (1999), the sizes are 
from Lawrence et al. (1986) . 



3.1.3. Radio lobe distance ratio Q 

Carilh et al. (1997) and Pentericci et al. (2000) have ob- 
tained 4.7 GHz and 8.2 GHz VLA observations of 64 
HzRGs. For the objects for which they identified a radio 
core, we measured the radio lobe distance ratio Q from 



their 8.2 GHz contour plot^. We also measured the Q val- 
ues in the sample of Best, Rottgering & Lehnert (1999), 
and took published values for 3CR galaxies from Best et 
al. (1995) and McCarthy et al. (1991). For small and faint 
sources this parameter is difficult to determine, which lim- 
its the possibility to examine the dependency of Q, espe- 
cially at the highest redshifts, where the radio size of the 
sources in our sample is smaller. 

3.1.4. Other radio parameters 

We considered to use the lobe flux density ratio 
R (e.g., McCarthy, van Breugel & Kapahi 1991a), which 



could provide a rough indication of orientation effects. 
However, too few published high resolution radio data are 
available at present to perform a statistically significant 
study of the correlations involving this parameter. 

3.2. Spectroscopic parameters 
3.2.1. Line fiuxes 

We used published line fiuxes throughout this paper (see 
Appendix A for references). The apertures used to extract 
the one-dimensional spectra were generally chosen to in- 
clude all the fiux of the most spatially extended emission 
line. The spread in aperture sizes and depths of the spectra 
will give rise to an increased scatter when comparing line 
luminosities from different samples. Some objects were ob- 
served under non-photometric conditions, and so we only 
include these sources for line-ratio studies. We have also 
determined 5a upper limits to the flux of N V in 45 ob- 
jects. Because line ratios are least sensitive to uncertainties 
in the relative flux calibration, we shall emphasise studies 
involving line ratios at the expense of studies involving 



^ Q is defined by McCarthy, van Breugel & Kapahi (1991) 
as the ratio of the distance from the radio core to the more 
distant radio hotspot divided by the distance from the radio 
core to the closer of the radio hotspot. 
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detailed kinematics of the emission lines, which can only 
be determined from the brightest line, i.e., Lya. 

3.2.2. Lya asymmetry 

The Lya emission in HzRGs often shows absorp- 
tion profiles caused by H i surrounding the radio 



galaxy, (e.g., van Ojik et al. 1997, Rottgering et al. 1995, 
Dey 1999| , |De Breuck et al. 1999[ ). This absorption pref- 
erentially occurs on the blue side of the emission line, 
which leads to a characteristic triangular shape of the Lya 
line in two dimensional spectra. Such asymmetries in Lya 
are often also observed in other objects with strong Lya 



emission at very high redshift (e.g., Dey et al. 1998), and 
might also have a contribution from intervening H i ab- 
sorbers along the cosmological line of sight. 

While high resolution spectroscopy of these Lya lines 
allows one to determine characteristics of the Hi ab- 



sorber ( [Rottgering et al. 1995| , [van Ojik et al. 1997| ), this 
is rarely possible with their discovery spectra, which have 
typical resolutions of a few hundred. To derive some in- 
formation on the H i absorption from these low-resolution 
spectra, we therefore introduce a parameter which mea- 
sures the relative flux on the blue and red sides of 
the assumed systemic redshift. Because absorption near 
the maximum of the Lya emission can shift the ob- 
served peak of the profile by as much as 500 km s~^ 
(Rottgering et al. 1997), we have to use other emission 
lines to determine the systemic redshift. The only other 
strong lines in HzRGs observed simultaneously with Lya 
are C IV, He H, and G HI]. From these. He II is the most 
appropriate line to use, because unlike C IV, He II is not 
a resonant line, so the profile should be less affected by 
absorption. 

We obtained one-dimensional spectra from the sam- 
ples of Rottgering et al. (1997), Stern et al. (1998) , and 
De Breuck et al. (2000b) . In 31 objects, we were able to 
obtain a good Gaussian fit of the He II line. We approxi- 
mated the error in the position of the peak of the Lya line 
as the quadratic sum of the error in the He II Gaussian 
fit and the error in the wavelength calibration, the latter 
being approximated as a quarter of the dispersion of ap- 
propriate spectrum (the resolution of the different spectra 
varies from ~10A to ~25A). 

We define the Lya asymmetry parameter 



A 



Fbi-u 



where Fuue is the flux within 2000 km s^^ blue- ward of 
the systemic Lya wavelength and Fred is the flux within 
2000 km s~^ red-ward of the systemic Lya wavelength 
(see Fig. ^ for an example). Because in the lowest reso- 
lution spectra, these 2000 km s~^ intervals often include 
only a few dispersion elements (pixels) of the spectrogram, 
we include only the percentage of the flux in the lowest 
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-1.0 -0.5 0.0 0.5 
Lya asymmetry parameter Al^^ 
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Fig. 4. Distribution of the Lya asymmetry parameter 
A-hya- Note the preference for blue- ward absorption (neg- 
ative Ai^ya values). 



and highest bin that lies within the 2000 km s~^ inter- 
val. A value of ^Lya^^l means total absorption of the 
flux on the blue side, a completely symmetrical profile has 
^LyQ=0, and positive Ai^ya values indicate absorption on 
the red side. We calculated the range of allowed ^Lya val- 
ues from the error in the predicted systemic wavelength 
of Lya. Because the peak sometimes falls on a steep side 
of the observed Lya profile, this can lead to asymmetric 
errors in Aj^ya ■ We excluded three sources where the range 
in ALyawas larger than 1.5 due to their large uncertain- 
ties. We could accurately determine the Aj^ya parameter 
for 25 HzRGs. The distribution (Fig. ^) shows a strong 
peak around Ai^ya=0 (no asymmetry) and a secondary 
peak around ALya=— 0.6 (blue absorption), but no sources 
with strong red-ward absorption. This is consistent with 
the trend seen from high resolution spectroscopy, where 
the median velocity of the Hi absorbers is blue-shifted 
by 100 km s~^ with respect to the peak emission redshift 



( |van Ojik et al. 1997] ). 

For 12 HzRGs, the resolution and signal to noise of 
the Lya profile was sufficient to extrapolate over the ab- 
sorption profiles while fitting the emission with a single 



Gaussian profile (see De Breuck et al. 1999 for an exam- 
ple). In cases where the profile is clearly non-Gaussian, we 
also included a Voigt function. We again used the He II 
redshift to predict the peak position of the Gaussian Lya 
emission (in two cases a variation of 1-2A provided a bet- 
ter fit, which is consistent with the wavelength calibra- 
tion errors). We varied the peak flux and FWHM of the 
Gaussian until we obtained a good flt on one side of the 
proflle. We use these crude Gaussian models to determine 
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4500 4550 4600 6120 6140 6160 6180 6200 
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Fig. 3. Example of the determination of the Lya asymmetry parameter Ai^ya. The right panel shows the position of 
the He II line (dash-dotted line) determined by a Gaussian fit with the uncertainties indicated by dotted lines. The 
left panel shows the position of the peak of Lya as predicted from the He II fit, with the errors indicated. The left and 
right hatched areas indicate the 2000 km/s intervals used to calculate the flux on both sides of the assumed systemic 
redshift. The ^Lya derived from this spectrum is ^Lya^ —0.49, with the range of uncertainties going from —0.58 to 
-0.39. 



how much flux was emitted before the H i absorption (inte- 
grating out to three times the FWHM of the Gaussian fit). 
We used this total flux to calculate an approximate frac- 
tion of the Lya flux that is absorbed. Figure ^ compares 
this fraction with ^Lya- As expected, we find a strong 
correlation indicating that (i) our Ai^ya parameter is a 
good approximation of the absorbed flux and (ii) that 
even in the most asymmetric profiles, the observed Lya 
flux is only diminished by ^ 50%, although we cannot ex- 
clude the possibility that sources with higher absorption 
fractions are missing from the samples because they are 
too weak to be detected. A limitation of our method is 
that it is insensitive to absorption which is centered near 
the peak of the emission, but such cases are rather rare 
( |van Ojik et al. 1997| ). 

For two of the highest redshift HzRGs where He II 
falls outside the optical window, we use the fit in Fig- 
ure H to calculate the approximate ^Lya value. This is 
particularly important if we want to study the redshift 
evolution of ^Lya- As a final note, it should be stressed 
that this procedure provides only a rough measure of the 
Lya kinematics. The main goal of the ^Lya parameter is 
to examine the existence of strong trends with redshift, 
or parameters of the radio source. A detailed study of the 
kinematics is impossible with these low-resolution data, 
and is beyond the scope of this paper. 



4. Correlations between parameters 




-0.6 -0.4 -0.2 0.2 

Lya asymmetry parameter A^y^, 

Fig. 5. The fraction of Lya flux absorbed plotted against 
the Lya asymmetry parameter (see text for details). Note 
than no more than 65% gets absorbed, and that the asym- 
metry parameter is a good measure of the total amount 
of absorption, as indicated by the linear fit. 

4-1. Selection effects 

The main problems we encounter when searching for cor- 
relations between the different parameters in our sample 
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of HzRGs are the selection effects resulting from the search 
techniques used to find them and technical limitations of 
the spectrographs. A well known example of the first ef- 
fect is the Malmquist bias in flux density limited samples. 
As discussed in §2, the addition of several filtered surveys 
with lower flux density limits alleviated this bias, espe- 
cially at 2 < z < 4, the range we shall concentrate on 
(Fig. [I]). However, comparisons with z < 2 radio galaxies 
will still be strongly affected by Malmquist bias. 

More subtle selection effects are due to the wavelength- 
dependent sensitivity of the optical spectrographs. Weak 
emission line objects could well have been escaped de- 
tection at certain redshifts where the brightest lines 
are not observable. For example, fewer sources are 
known in the 1.5 < z < 2 "redshift desert" where 
[O n] A 3727 has shifted out of the optical window and 
Lya has not yet entered, and the redshift will have to 
be based on weaker lines like C IV, He H C HI] or C II] 
(e.g., ^tern et al. 1999| ). At very high redshift (z > 3), 
weaker lines red-ward of Lya will shift to wavelengths 
where measurement is difficult due to fringing of the CCD 
and bright OH sky-lines (e.g., Osterbrock fc Martel 199^ . 



lent widths with survival analysis methods, we shall use 



the ASURV software package ([sobe & Feigelson 1990, 
Lavalley, Isobe & Feigelson 1992), which is incorporated 



This should have only a limited effect on the line ratio di- 
agnostics, because the spectra of z > 3 radio galaxies are 
generally of higher quality than those at 2 < z < 3: the 
highest redshift spectra were usually obtained with larger 
aperture telescopes, or the integration time was extended 
to search for confirming lines. 

4-2. Statistical tests 

Most previous studies have each concentrated on spe- 
cific correlations between a few specified parameters in 
HzRGs. Here, we shall examine correlations between all 
of our measured parameters. Some apparent correlations 
between two parameters might be due to the correlation 
of one or both of these parameters with a third parameter 
(either by selection effects or by a real correlation). As ar- 
gued by Macklin (1982), the statistical test that de-couples 
such dependent correlations, while retaining as much in- 
formation as possible, is the Spearman partial rank cor- 
relation coefficient. We shall apply this test to determine 
if some correlations are not due to selection effects in one 
or both of the parameters. 

In the first approximation, we simply exclude lower 
or upper limits from this analysis. For most parameters, 
this will not have a serious effect on the results, be- 
cause less than 5% of the data are lower or upper lim- 
its. However, 12% of the linear sizes are upper limits, and 
an even higher percentage of the equivalent widths are 
lower limits (from 11% for W^lfl to 35% for Wl^'J). To 
treat these limits, or "censored data", in a statistically 
meaningful manner with minimal loss of information, a 
class of statistical me thods called "survival analysis" has 
been developed (e.g., [sobe, Feigelson fc Nelson 1986 ). To 
test the correlations involving radio size and cquiva- 



into the NOAO reduction package IRAF. Three different 
correlation tests within ASURV are of importance to this 
work: Cox's hazard model, generalized Kendall's r, and 
generalized Spearman's p. We shall take the variables as 
correlated if all three of these models give significance lev- 
els^ >95%. One of the deficiencies of survival analysis for 
astronomical purposes is that upper limits are assumed 
to be precisely known. Especially the lower limits to the 
equivalent widths are very uncertain due to the extreme 
faintness of the continuum emission. Survival analysis is 
not as effective in handling dependent correlations in the 
way the Spearman partial rank analysis can. We shall 
therefore only use survival analysis to check whether the 
addition of censored data changes the significance of a 
correlation which was found from our previous analysis. 



4-3. Spearman rank analysis 

To investigate every possible dependent correlation, we 
first calculated the Spearman rank correlation coefficient 
and the associated significance levels for all 153 possible 
combinations of source parameters in our HzRG sample. 
The results are presented in appendix B. The correlation 
coefficients and number of available parameters are tabu- 
lated in the lower left half of the Table, and the significance 
levels in the upper right half. Wc also calculated these val- 
ues for a. Hq ~ 65 km s~^Mpc^^, qo = 0.15 cosmology, but 
this changed the values of the correlation coefficients by 
less than than 0.1. 

We find 33 correlations with significance levels >99%, 
and 10 additional which are between 95% and 99% signif- 
icant. These correlations are: 

• Redshift z with radio size D, radio power P325, radio 
spectral index a and line luminosity Lciv, -^HcII or £ciii]- 

• Radio spectral index a with radio size D. 

• Radio spectral index a with radio power ^325- 

• Radio spectral index a with line luminosity (increasingly 
stronger correlation for weaker emission lines). 

• Radio spectral index a with equivalent width of C IV, 
He II and C HI]. 

• Radio size D with Lya asymmetry ^Lya, and £hcII- 

• Radio lobe distance ratio Q with equivalent width W[y^^ 
and W^lll. 

• Radio power P325 and P1400 with all four line luminosi- 
ties. 

• Line luminosities of Lya, C IV, He II and C HI] with 
each other. 



Note that for consistency with previous statistical work on 
HzRGs, by significance level we mean the significance level of 
the correlation, and not the probability of the null hypothesis 
(no correlation) being true, which is generally used in survival 
analysis. 
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• Lya luminosity with both Lya equivalent width W£y^ 
and Lya line width AwLya- 

• Equivalent widths of Lya, C IV, He II and C III] with 
each other. 

• Equivalent widths of fainter lines (He II and C HI] ) with 
line widths of Lya and to a lesser extent fainter lines. 

We shall now examine if some of these correlation are 
due to selection effects, dependent correlations with other 
parameters, or are removed if sources with upper/lower 
limits are included. 

4-.4- Correlations influenced or caused by selection effects 

4.4.1. Radio size and power vs. redshift 

The dependence of radio size D on redshift z 
and radio power P325 has been examined by 
a number of authors (e.g ., Neeser et al. 1995 , 
Blundell, Rawlings fc Willot 19991) . The range of of- 
ten contradictory results is due to several selection effects 
(e.g., in selection frequency) in the samples used to 
examine this correlation. Our sample of radio galaxies 
is too inhomogeneous and incomplete to address this 
question. For example the Malmquist bias in radio power 
(Fig. |l|) will dominate the P325 — z correlation. For z > 2, 
the probability that z and P325 are uncorrelated drops 
from 1.15 X lO^'^^ to 1.6 x 10^"*, indicating the Malmquist 
bias is strongly decreased, but not removed (see Fig. 
|l|). While part of the correlation might be due to an 
underlying real correlation (e.g., between z and D), we 
consider the effect of the selection effects too large and 
will not examine the significance of these relations. 



4.4.2. Radio spectral index vs. radio size and radio power 

In some samples designed to find HzRGs, the radio size 
D has been used as an additional "high redshift filter" in 
combination with a spectral index cutoff. This will lead 
to selection effects in the redshift-spectral index relation, 
masking a possible real correlation. 

The relation between spectral index and radio power 
is not independent, because we used the spectral index to 
calculate the rest-frame radio power of our sources. 



4.4.3. Radio spectral index vs. line luminosity and 
equivalent width 

The luminosities of the UV emission lines and rest-frame 
equivalent widths appear to be correlated, with the weak- 
est line (C HI]) showing the strongest correlation. How- 
ever, due to the way we constructed our sample of radio 
galaxies, the sources with the flattest spectral indices are 
all at relatively low redshift (z < 2). Because the line lu- 
minosities are strongly correlated with redshift (see next 
section), we only find the low-luminosity lines in the low- 
redshift flattest spectrum sources of our sample. 
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Fig. 6. Rest-frame equivalent width plotted against radio 
spectral index. Note that the trend for redder lines (He II 
and C III]) is more offset towards flatter spectral indices 
than for the blue lines (Lya and C IV). This arises because 
at higher redshift (i.e., steeper spectral index) the bluest 
lines either shift into less sensitive parts of the observable 
window, where they can only be detected if they have high 
equivalent widths. 



The apparent correlation between radio spectral index 
and rest-frame equivalent width is also influenced by sub- 
tle redshift selection effects: the increasingly bluer lines 
will only be detected at higher redshifts due to the fixed 
wavelength range of the optical spectrographs. At the low- 
est redshifts, these lines will be in the bluest, least sen- 
sitive part of the CCD, and they will only be detected 
when they have a high equivalent width. This effect is 
obvious in Fig. |^: the apparent correlation between rest- 
frame equivalent width and spectral index seems to shift 
towards flatter spectral indices in the redder lines, which 
are only observed in the lower redshift sources which have 
flatter radio spectral indices. 

4.4.4. Emission line luminosity vs. redshift 

The apparent correlation between line luminosity and red- 
shift z is a less known selection effect. The correlation 
is weak or absent for Lya, but increasingly stronger for 
C IV, He II, and C HI]. The reason for this is obvious 
in Figures ^ and ^, where we plot the emission line lu- 
minosities against redshift together with a curve which 
denotes the luminosity corresponding to an emission line 
with a flux of 1.0 x 10~"'^^erg s~^ cm~^A~-'^, which is close 
to the minimum detectable level of most of the obser- 
vations. We clearly find a line-luminosity redshift degen- 
eracy, as found for the radio powers. This degeneracy 
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Fig. 7. Emission line luminosity for C IV, He II, C III] 
and C II] plotted against redshift. The plotted line 
indicates the luminosity as a function of redshift at 
which a line with an emission line flux of 1.0 x 
10~"'^^erg cm~^A~^, near the detectability of 3-4m 
class telescopes, would be detected. Note that the higher 
wavelength lines are increasingly weaker, and lie closer to 
the detection limit. The lack of sources below this limit 
leads to an artificial redshift dependence of the line lumi- 
nosity. 

Table 2. Results of survival analysis on the correlations 
between line luminosity and equivalent width. The signif- 
icance indicates the probability a correlation is detected. 



Vai 


lablc 


Percentage 




Signitican 




Independent 


Dependent 


Censored 


Cox 


Kendall 


Spearman 






36% 


34.71% 


85.00% 


81.36% 


^CIV 




18% 


99.92% 


99.78% 


99.65% 


^Hell 


yi^rest 
HcII 


11% 


99.51% 


99.25% 


99.31% 


^CIII] 


^rest 
^CIIIl 


10% 


70.19% 


98.29% 


97.85% 



becomes more pronounced as the lines become weaker. 
The Lya line is generally an order of magnitude brighter 
than the other UV lines. In addition, the two highest red- 
shift radio galaxies have very weak Lya (TN J0924— 2201 
at z = 5.19 and VLA J 123642-H621331 at z ^ 

Waddington et al. 1999D, 



4.42; van Breugel et al. 1999 



and have been detected at flux levels an order of magni- 
tude fainter than what has been attempted for the other 
lines. The lines red-ward of Lya are increasingly fainter, 
and lie much closer to the detection limit, and will there- 
fore more easily be missed. 

4.4.5. Emission line luminosity vs. equivalent width 

The correlations between the equivalent width and line lu- 
minosities are also dominated by selection effects, namely 



Fig. 8. Emission line luminosity for Lya plotted against 
redshift. The line indicates a line with an emission line flux 
of 1.0 X 10~^^erg cm~^A~^, near the detectability of 3- 
4m class telescopes. Note that there is no artificial redshift 
dependence, as for the lines in Figure ^ 



the difficulty to detect objects with weak emission line 
fluxes and high equivalent widths. These sources probably 
do exist, as the objects for which there are lower limits to 
the equivalent widths frequently have low line luminosi- 
ties. When we apply survival analysis models to the com- 
plete dataset (Table ||), we find that the correlations with 
Lya and C III] are no longer significant. Given that both 
the detections and lower limits to the equivalent widths in 
these weak lines are highly uncertain, we interpret these 
correlations as too uncertain to be trustworthy. 

4.4.6. Lya line width vs. equivalent width of He II and 
C III] 

A final correlation which we consider dominated by se- 
lection effects is the correlation of equivalent widths of 
fainter lines (He II and C III]) with line widths of Lya. 
When the two sources with AuLya > 3000 km s~^ are not 
considered, the significance of these correlations drops to 
< 95%. Such sources probably have a broad component 
detected in Lya, but not in the weaker lines. 



4-5. Possible correlations 

A few parameters, for which limited data exist appear to 
show correlations. One set of such correlations are those 
between the radio lobe distance ratio Q and rest-frame 
equivalent widths of the UV lines. From the plots of these 
four correlations (Fig. ^), we find that there is a dearth 
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Fig. 9. Equivalent widths of the strong UV hnes plotted 
against radio lobe distance ratio Q. 



of highly asymmetric radio sources with large equivalent 
widths. 

One possible explanation of this effect can be found in 
the observing techniques used to determine the redshifts of 
HzRGs. In some objects where the correct identification of 
the host galaxy is uncertain, the spectroscopic slit will be 
positioned such as to increase the chances to detect strong 
line emission which can be used to determine the redshift. 
Most often, this means that the slit is aligned with the 
radio emission. In relatively symmetrical sources with a 
slight bending of the radio lobes, this could mean that the 
slit will not be perfectly centered on the host galaxy. In 
highly asymmetric sources, the slit will be more often cen- 
tered near the presumed optical identification, which will 
increase the chance that host galaxy continuum emission 
is included, and lead to a lower equivalent width. Based 
on only ~20 sources, we therefore consider this correla- 
tion suggestive. A sample of radio galaxies observed in a 
consistent way would be needed to fully examine the sig- 
nificance of this correlation. 



Variables 


Number 


rpL 




a 




67 


0.35 (0.38) 


0.37 (0.38) 


3.06 (3.14) 


^325 1 i'Civ 


64 


0.40 (0.28) 


0.29 (0.28) 


2.40 (2.27) 


-P325 ) .^Hell 


61 


0.35 (0.24) 


0.25 (0.24) 


1.94 (1.91) 


-P325 . -f-CIII] 


79 


0.47 (0.30) 


0.31 (0.30) 


2.78 (2.68) 


^"325 1 -'^CII] 


42 


0.68 (0.73) 


0.70 (0.70) 


5.29 (5.39) 


^^325 1 ^Mgll 


30 


0.62 (0.65) 


0.56 (0.56) 


3.25 (3.26) 



Table 3. Spearman partial rank correlation analysis for 
the correlations between radio power P325 , line luminosity 
of the bright UV lines and redshift. rp^ is the Spearman 
rank coefficient of the radio power and line luminosity, and 
TpL.z is the Spearman partial rank coefficient, taking the 
selection effects with redshift into account, a is the signif- 
icance of the partial rank correlation, which is equivalent 
to the deviation from a unit variance normal distribution 
if no correlation is present. Values in brackets are for a 
Hq = 65 km s^^Mpc^^, qo^O.lS cosmology. 



4.6.2. Radio power vs. emission line luminosity 

To examine the correlation between radio power and emis- 
sion line luminosity, we use the Spearman partial rank 
coefficient, including redshift as the third parameter, and 
thus taking the Malmquist bias of both parameters into 
account. Table |3| presents the results of this analysis. Fig- 
ure ^ shows graphical representations of the correlations. 

Table ^ shows that the selection effects with redshift 
have only a minor influence on the correlation between the 
radio powers and Lya C II] and Mg II luminosity, but do 
artiffcially strengthen the correlations with C IV, He II, 
and C HI]. This is because the radio power - redshift de- 
generacy becomes less pronounced, or absent at z > 2, 
where the UV lines shift into the optical spectroscopy win- 
dow, and because of the increasingly more important red- 
shift selection effects with weaker emission lines, as dis- 
cussed in §4.3. We have examined this correlation using 
the radio power at 325 MHz and 1.4 GHz, and found only 
minir differences, with the P325 yielding slightly stronger 
correlations. We prefer to use the P325 because the higher 
frequency might have a larger contribution from a Doppler 
boosted core (see e.g., Blundell, Rawlings & Willot 1999). 

We find there is a > 2cr probability that the corre- 
lations between radio power and line luminosity do not 
arise from the correlations of individual parameters with 
redshift. We shall discuss this further in §6.2. 



4-. 6. Probable correlations 

4.6.1. Line luminosities and equivalent widths 

We find that the line luminosities of the four bright UV 
lines are all strongly correlated with each other. We also 
find a similar result for the rest-frame equivalent widths. 
This indicates that the different emission lines in HzRGs 
are powered by the same mechanism. We shall discuss this 
further in §5. 



4.6.3. Lyck asymmetry vs. Radio size, and vs. redshift 

The Spearman rank analysis (appendix B) indicates that 
there is an anti-correlation between radio size and Lya 
asymmetry Ai^ya with a confidence level of 99.86%. We 
perform a Spearman partial rank analysis, because of the 
possibility of this arising artificially from the correlation 
between redshift and radio size D. 

From this analysis (Table |^), we find that Ai^ya ap- 
pears to be correlated with radio size and redshift inde- 
pendently. From Figure |ll|, we see that all seven HzRGs 
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Fig. 10. Line luminosity for Lya, C IV, He II, C III] 
C II] and Mg II plotted against radio power at 325 MHz. 
Crosses represent sources at z < 3, triangles those at 
z > 3. 
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Fig. 11. Lya asymmetry parameter Alyq plotted against 
radio size. Error bars reflect the uncertainty in the wave- 
length calibration and fitting of the He II line, and need 
not be symmetrical. Note that strong blue Lya absorption 
only occurs in smaller radio sources. 



dicated, which all show strong blueward absorption. His 
results are not only consistent with the ALy^— D trend, 
but also show stronger absorption for the z > 4 radio 
galaxies, as is also seen in TN J1338— 1942 at z = 4.11 
(De Breuck et al. 199E). We discuss this result further in 
§6.4. 



X; V 


Number 














25 


-0.60(- 


-0.59) 


--0.64(- 


-0.64) 


3.48 (3.44) 




25 


0.33( 


0.33) 


0.41( 


0.09) 


2.02 (2.06) 



Table 4. Spearman partial rank correlation analysis for 
the correlations between radio size _D, Lya asymmetry 
^Lyaj and redshift z. r^y is the Spearman rank coeffi- 
cient, and Txy^z is the the Spearman partial rank coeffi- 
cient, taking the selection effects with the third parame- 
ter into account, a is the significance of the partial rank 
correlation, as in Table Values in brackets are for a 
Hq = 65 km s^^Mpc^^, go=0-15 cosmology. 



4.7. Summary of correlation analysis 

After consideration of all selection effects, we find 22 real 
correlations, viz. 

• Radio size D and redshift z with Lya asymmetry Ai^ya ■ 

• UV line luminosity with radio power (8 correlations). 

• Line luminosities of Lya, C IV, He II and C HI] with 
each other (6 correlations). 

• Equivalent widths of Lya, C IV, He II and C HI] with 
each other (6 correlations). 



with Ai^ya> 0.3 have radio sizes Z? < 60 kpc. This con- 
firms the results of van Ojik et al. (1997), who found from 
high resolution spectroscopy that H i absorption preferen- 
tially occurs in radio galaxies smaller than 50 kpc. 

The correlation with redshift is less significant, but wc 
lack sufficient data at z > 3, because our method requires 
the detection of the He II line, which starts shifting out of 
the optical band at these redshifts. However, Dey (1999) 
shows Lya velocity profiles of four z > 3.5 radio galax- 
ies with the redshift predicted from He II or C HI] in- 



5. Emission line ratios 

The relative intensities of emission lines provide a powerful 
tool to examine the physical state of the extended emis- 
sion line gas in HzRGs. The (rest- frame) UV emission lines 
from different elements (mainly hydrogen, carbon, helium, 
and nitrogen) and the three ionization states of the car- 
bon lines (C II], C III], and C IV) can provide valuable 
information on the relative abundance of these elements, 
and can distinguish between ionization mechanisms. 

An efficient way to compare the data with model pre- 
dictions is to use line-ratio diagrams. For HzRGs, the 
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rest-frame optical lines which are normally used for such 
studies have shifted to the near-IR, and we have to use 
the rest-frame UV-lines that have shifted into the op- 
tical. VM97 and ADT98 have used line ratios involving 
the C IV, He II, C III] and C II] lines to compare HzRG 
spectra (from Rottgering 1997) with the predictions from 
photo-ionization and shock models. Our HzRG sample 
contains almost three times as many sources as the sam- 
ples in VM97 and ADT98: 58 sources have simultane- 
ously observed C IV, He II, and C HI] lines, with redshifts 
1.2 < z < 3.8, radio sizes D < 365 kpc, and radio powers 
34.5 < log(P325) < 36.5. In this section, we shall examine 
the dependence of the line ratios on these parameters, and 
on the physical parameters determined from the shock and 
photo-ionization models. Because we have seen in §4.4 that 
all of these parameters have an artificial redshift depen- 
dence, we shall first disentangle the mutual dependences 
of the line ratios on these parameters. After introducing 
the theoretical models, we shall compare the data with 
the range of model predictions, and examine the ratios in- 
volving Lya and two rest-frame optical lines determined 
from near-IR spectroscopy. 

5.1. HzRG UV line ratios 

In §4.4, we found an increasing dependence of the C IV, 
He II, C HI] and C II] line luminosities with redshift. This 
selection effect will also affect the line ratios. To properly 
examine the trends involving line ratios, we used a four- 
parameter Spearman partial rank analysis (Table |^) . 

We find that only one of the correlations with more 
than 20 elements has more than 2a significance, indicat- 
ing that redshift, radio size or power do not have a strong 
influence on most emission line ratios. The only one of 
these correlations that is significant is the one between the 
C II]/C HI] ratio and radio size D, which BRLOO reported 
as evidence that small sources in their sample of 14 3CR 
galaxies at z '^l are ionized by shocks and larger sources 
by photo-ionization. Our sample, which includes all their 
sources, is twice as large, but only adds relatively small 
sources with D < 150 kpc. The Spearman partial rank 
analysis shows that the highest probability of correlation 
of C II]/C HI] is indeed with D, and not with z or P325. 
The results of BRLOO were not strongly affected by selec- 
tion effects, because their 3CR subsample has a limited 
redshift range 0.7 < z < 1.2. The addition of six upper 
limits to the radio size (Fig. [ij) is consistent with the 
C H]/C HI] - D correlation because no sources > 50 kpc 
with C II] stronger than C III] are found. 

Stern et al. (1999) (hereafter S99) reported a possi- 
ble correlation between the C IV/C HI] ratio and radio 
power. We find no evidence for a correlation of C IV/C III] 
with P325, eventhough S99 used 25 of our 43 data points. 
This is partly due to the exclusion of five sources (of which 
three are from S99) with only upper limits to the ra- 
dio size. These three excluded sources are the only ones 
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Fig. 12. The C II]/C HI] ratio plotted versus radio size. 
Note the absence of large sources with strong C II] . 

where C HI] is stronger than C IV. If we retain these five 
sources and replace their upper limits to the radio size by 
detections^, we find a Spearman partial rank coefficient 
frP.zD = 0.22 and significance level Drp^zD = 1.35. Fur- 
ther differences are caused by our different determination 
of the radio power (§3.1.1). It is clear that a larger and 
better defined sample with multi-frequency radio data to 
determine the radio power and deep spectroscopy to break 
the artificial redshift dependence of the line luminosities 
is needed to examine whether the C IV/ C HI] line ratio 
is correlated with redshift, radio size, or radio power. 

To summarize, the only strong correlation between the 
UV line ratios and radio properties is a dependence of the 
C II] /C HI] ratio on radio size D. We find no evidence 
that the radio power and ionization state are correlated. 

5.2. Shock and photo-ionization models 

Of great importance in the study of the nature of the ex- 
tended emission line regions in HzRGs is the mechanism 
responsible for producing the observed ionization state. 
The most likely mechanisms include photo-ionization from 
a central AGN, and shock ionization from the radio jets 
propagating through the interstellar medium. In this sec- 
tion, we describe three theoretical models that can predict 
the line ratios for these processes. 

To calculate the predictions from shock ionization, 
we use the models of Dopita & Sutherland (1996; here- 
after DS96). Their models consist of two components: 
(i) a down-stream shock component, dominated by the 
radiative cooling of the gas behind the shock, and (ii) 

* As discussed in §4.2, we cannot use survival analysis on 
dependent variables. 



Carlos De Breuck et al.: Statistical study of emission lines from HzRGs 



13 



Line ratio 


Number 




z 






D 






P325 






Trz.DP 


a 


TrD 




a 


TrP 


rrP,zD 


a 


Lya / C IV 


41 


-0.13 


-0.06 


-0.34 


0.11 


0.08 


0.47 


-0.09 


-0.05 


-0.33 


LyQ / He II 


32 


-0.24 


-0.06 


-0.30 


0.15 


0.16 


0.83 


-0.35 


-0.32 


-1.71 


Lya / C HI] 


27 


-0.19 


-0.20 


-0.93 


-0.16 


-0.23 


-1.10 


-0.19 


-0.17 


-0.81 


C IV / He II 


48 


-0.10 


-0.13 


-0.87 


-0.05 


-0.08 


-0.53 


0.02 


0.08 


0.51 


C IV / C III] 


43 


0.08 


0.09 


0.55 


-0.01 


-0.00 


-0.01 


-0.05 


-0.06 


-0.40 


C IV / C II] 


11 


0.19 


0.15 


0.38 


0.28 


0.53 


1.45 


-0.74 


-0.80 


-2.71 


C III] / He II 


45 


-0.24 


-0.25 


-1.65 


0.03 


-0.03 


-0.18 


0.03 


0.09 


0.57 


C II] / He II 


14 


0.32 


0.19 


0.57 


-0.63 


-0.65 


-2.31 


0.31 


0.58 


1.97 


C II] / C HI] 


29 


-0.41 


0.08 


0.37 


0.69 


0.60 


3.39 


0.03 


-0.08 


-0.39 



Table 5. Spearman partial rank correlation analysis for the correlations between nine UV line ratios and redshift z, 
radio size and radio power P325. Vrz is the Spearman rank correlation coefficient between the line ratio r of the 
lines in the first column and z, while rrz,DP is the Spearman partial rank coefficient of that correlation, taking the 
possible correlations with D and P325 into account, a is the significance of the partial rank correlation, as in Table 2. 
The last six columns similarly examine the correlations with radio size and radio power. 



an up-stream precursor component, dominated by the 
photo-ionization by the radiation from the shocked gas. 
Their models assume solar metallicities, and present a 
grid of model predictions for shock velocities in the range 
150 < V < 500 km s^^, and for values of the magnetic 
parameter < Bj^fn < 4/iG cm^/^, with B the pre- 
shock transverse magnetic field and n the pre-shock num- 
ber density. This magnetic parameter controls the effective 
ionization parameter in the down-stream emitting com- 
ponent of the shock since at high shock velocities, the 
transverse magnetic field limits the compression caused 
by the shock through a balance of the magnetic pres- 
sure of the cloud and the ram pressure of the shock. Fu- 
ture versions of these models extend the shock velocities 
out t o 1000 km s~^, and will include non-solar metallici- 
ties ( Bickncll ct al. 200C| ; Sutherland, private communica- 
tion). To derive the He II A 1640 fluxes (not included in the 
tables of DS96) from the He II A 4686 fluxes that are pro - 
vided, we assume a ratio of 9 (e.g., MacAlpine ct al. 1985| ). 
Note that the shock+precursor models of ADT98 do not 
include this factor for the precursor component so that 
their He II A 1640 values are underestimated by a factors 
of ~ 1 — 6 (Allen, private communication), whilst VM97 
plotted the shock and precursor models separately, which 
is inappropriate for HzRG spectra, since these are spa- 
tially integrated over both the shock and precursor gas. 

For the photo-ionization models, VM97 use the MAP- 
PINGS I code developed by Binette, while ADT98 use the 
MAPPINGS II code developed by Sutherland. The diflfer- 
ences between the two model predictions are minor for 
the relevant range of parameters. We use CLOUDY ver- 



sion C94 (Ferland 1996), and also find very similar results. 



providing another independent consistency check with the 
results of VM97 and ADT98. We calculate the model spec- 
tra for the same range of parameters as ADT98 to facil- 
itate comparison between the predictions of the different 
codes. We assume solar metallicities, and an ionizing con- 



tinuum which is a power law spectrum ((/)^ cx v") with 
lower and higher energy cutoffs of 0.01 Ryd and 100 Ryd 
(1.36 keV), but the models do not depend strongly on the 
exact values of these limits. We calculate spectra for a 
power law spectral index a — —1.5 which well matches 
the low redshift Seyfert spectra (e.g., Evans et al. 1999), 
and for a = -1.0, which VM97 find to match the HzRG 
data better. For the hydrogen density, we used a value of 
n — 100, as co mmonly used for the ex tended emission line 
regions (e.g., McCarthy et al. 1990 ) and a high density 
value of n = 1000 cm~"^. We vary the ionization parame- 
tei0 U from 0.001 to 0.1 in steps of 0.5 dex to cover the 
entire range of observed line ratios. 

Binette, Wilson & Storchi-Bcrgmann (1996, hereafter 
B96) present an alternative photo-ionization sequence by 
considering emission from two distinct cloud populations. 
In their model, the light from the photo-ionizing source 
first passes through a population of matter-bounded (op- 
tically thin) clouds, and subsequently strikes ionization- 
bounded (optically thick) clouds located further outwards. 
They produce a sequence by varying the the parameter 
^M/i; defined as the ratio of the solid angle subtended by 
the matter-bounded clouds to the solid angle subtended by 
the ionization-bounded clouds. They keep the power law 
spectral index of the ionizing continuum incident at the 
matter-bounded clouds constant at a = —1.3 and the ion- 
ization parameter constant at Umb = 0.04. By adding this 
ionization-bounded component, B96 solve several short- 
comings of the classical photo-ionization models, like the 
inability to produce strong high excitation UV lines for 
reasonable values of U inferred from the optical lines. Be- 
cause the ^M/i models require that the ionizing contin- 



^ Defined as U — (en) ^ {(j)„diy)/hi', with </!>„ the 
monochromatic ionizing energy fiux impinging on the cloud, vq 
the ionizing potential of hydrogen, n the total gas density at 
the front face of the cloud, c the speed of light, and h Planck's 
constant. 
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uum incident on the ionization-bounded clouds first has 
to be "filtered" by the matter-bounded clouds, the val- 
ues of Ajvi/i cannot strictly be less than unity. However, if 
the matter-bounded clouds are obscured along the line of 
sight, we can observe apparent Am/i<1 values. 



distinct regions in the galaxies. The fact that the ^m/i 
sequence provides the best fit in this diagram (see below) 
is consistent with this idea, as in this model we observe si- 
multaneously emission from clouds with different incident 
ionizing radiation. We return to this in §6.1. 



5.3. Comparison with the observations 

Figure ^ reproduces three diagnostic diagrams from 
VM97 and ADT98, showing the C IV, He H, and C III] 
lines, and a Carbon-only plot from ADT98, showing the 
ratios of C IV, C III], and C II]. 

5.3.1. Shock models 

From the left panels in Figure |l^, we find that the present 
shock or shock -|-precursor models cannot explain the ra- 
tios involving C IV, He II, and C III] , while roughly half of 
the HzRGs in the Carbon-only plot coincide with the high- 
est shock velocity (500 km s~^) models. In the plots with- 
out C II], the least discrepant predictions are the high- 
est shock velocity shock-|-precursor models. The disturbed 
kinematics of the emission line gas indicate that shocks 
with velocities >500 km s^^ are probably occurring in 



HzRGs (Bicknell et al. 2000), so it will be of great inter- 
est to compare the data with higher shock velocity models 
(Sutherland et al. , in preparation). Preliminary results 
indicate that the existing models cannot be simply ex- 
trapolated (Allen, private communication). Moreover, the 
DS96 shock models have only been calculated for solar 
metallicities, which are unlikely to be correct in the ex- 
tended emission lin e regions at the highest redshifts (§6.3, 
Binette et al. 200"o| ). Lower metallicities will lead to less 



efficient cooling and raise the temperature in the precur- 
sor area, leading to stronger emission. We conclude that 
the present shock models can only reproduce the C II] 
vs. C HI] ratio, but that extensions of the present models 
to higher shock velocities and lower metallicities will be 
needed to fully determine the importance of this ionization 
mechanism. 



5.3.2. Pure photo-ionization 

From the right panels in Figure |l^, we find that pure 
photo-ionization provides good fits to the C IV/He II and 
C IV/C HI] ratios, and a reasonable fit to the C HI] /He II 
ratio, but fails to explain the C II]/C HI] ratio. In the 
plots without C II] the majority of the data-points are 
bracketed by the models with power law spectral indices of 
the incident ionizing continuum of a = — 1 and a = —1.5. 
In the plot involving C II] the a = — 1 model provides the 
least discrepant predictions. All plots suggest ionization 
parameters — 2.5 < logio(U) < — 2 for the a — ~\ models 
and —2 < logio(U) < —1.5 for the a = —1.5 models. 

The poor fit in the Carbon-only plot indicates that 
the different ionization stages of Carbon originate from 



5.3.3. Am /I sequence 

The predictions for the A-^/i sequence in the models of 
B96 provide good fits to the data, but fails to reproduce 
the sometimes large spread of the points. This could be 
solved by changing the dust content, metallicity, spectral 
index or ionization parameter Umb of the ionizing contin- 
uum incident on the matter-bounded clouds, or by inter- 
nal reddening^. A value of Am/i~0.1 seems to best fit the 
data in all four diagrams, although the scatter ranges from 
^M/i~0.02 to ^M/i^lOO. A value of ^m/i~0.1 means that 
as much as 90% of the matter-bounded clouds has to be 
obscured along the line of sight. In low redshift Seyferts, 
B96 and Evans et al. (1999) found values of A^/i shghtly 
higher than unity provided the best fit to the data. If the 
Am /I sequence is the correct ionization mechanism model, 
the matter-bounded clouds would be much more obscured 
at higher redshift, or more obscured in radio galaxies than 
in Seyferts. 

A possible way to obscure these regions could be pref- 
erential obscuration by dust. Dust masses up to IO^Mq 



have been detected in HzRGs (e.g., Archibald et al. 2000). 
If these large amounts of dust are located near the cen- 
tral parts of HzRGs, they could lead to obscuration of 
either the matter-bounded or ionization-bounded clouds. 
There are six radio galaxies in our sample with solid 
submm detections which also have published C IV, He II 
and C HI] measurements (4C41.17, Chini & Kriigel 1994, 
Dunlop et al. 1994; MG J1019-H0534, Cimatti et al. 1998a; 
4C24.28, 4C28.58 and 4C48.48, Archibald et al. 2000; 
TN J0121+1320, Reuland et al. , in preparation). These 
six sources are not preferentially located in any part of the 
diagrams, and do not fall in the low Am/i region. We con- 
clude that dust obscuration of the matter-bounded clouds 
is an unlikely explanation of the very low Am/i values. 
A variation of the Am/i sequence by changing the ioniza- 
tion parameter or power law spectral index incident on 
the matter-bounded clouds would be required to explain 
the high obscuration of the matter-bounded clouds and 
the observed scatter with these models. 



5.3.4. Summary of comparison with the data 

To summarize, we find that using only the line-ratio dia- 
grams involving C IV, He II, and C HI], the pure photo- 



® Galactic reddening will be negligible because most HzRGs 
have been identified from samples avoiding the Galactic plane, 
and some of the spectra in the literature have been corrected 
for reddening. 
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Fig. 13. Line ratio diagnostic diagrams for C IV, He II, C III], and C II]. The left panels show the shock (shaded 
grid) and shock+precursor (unshaded grid) models. Shock velocity increases along the lines from 150 km s~^ at 
the top or right to 500 km at the end of the curves, which show four values of the magnetic parameter < 
B I yfn < 4 fiG cm'^/^. The sequence shown by open triangles represents a composite model of a pure shock model with 
B / ^/n = 4 fiG cm'^/^, v =400 km and a photo-ionization model with power-law spectral index a = —1.0, hydrogen 
density n = 100 cm~^ and ionization parameter Logio(?7) — —2.25. The triangles represent different fractions of each 
model in steps of 20% (see text for details). 

The right panels show the photo-ionization models. We show four photo-ionization sequences, with the values of 
logio(C/) every 0.5 dex. The top shaded grid represents two models with power law spectral index a — —1, the bottom 
grid models are for a = —1.5. The boundaries of the shaded sequences are models for a hydrogen density n = 100 (left 
or bottom) and n — 1000 (top or right). The single curve is the A^/i sequence of B96, with values 0.01 <Am/i< 100 
and tickmarks every 0.2 dex. 

The data are repeated in both panels. Open circles represent sources at z < 2, circled -I— signs at 2 < z < 3, and filled 
circles at z > 3. The asterisks represent the two nearby Seyfert galaxies observed by Evans et al. (1999). 



ionization models provide the best fit to the HzRG data, 
but these models clearly fail to reproduce the observed 
C II]/ C III] ratio, which can be well fit by high shock 
velocity models. The models that use a combination of 



matter-bounded and ionization-bounded clouds seem to 
provide a reasonable fit to all four lines, but the spread 
of the points does not seem to follow the Am/i sequence, 
and suggests other parameters dominate the intrinsic dif- 
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Fig. 13. , continued. 



ferences in the ionization levels of the individual HzRGs. 
However, the combination of at least two zones with dif- 
ferent ionization continua seems to be required to explain 
the emission line ratios in C IV, He H, C HI], and C II]. 

In the remainder of this section, we shall compare these 
model predictions with some less commonly used emission 
line ratios. In §6.1, we shall return to the multiple zone 
ionization mechanisms. 



5.4- Diagnostic diagrams with Lya 

In almost all HzRGs, the brightest UV line is Lya. This 
line has not been frequently used in diagnostic line-ratio 
diagrams, because it is highly sensitive to resonance scat- 
tering inside and around the excited clouds, and to ab- 
sorption by dust. 



Villar-Martm, Binette & Fosbury (1996) have exam- 
ined the effects of resonance scattering and dust on UV 
lines in HzRG. They found that cases where the Lya emis- 
sion is weak relative to the other UV lines can better be ex- 
plained by geometrical (viewing angle) effects rather than 
by large amounts of dust, which would also affect the other 
UV lines, notably C IV. Another factor affecting the total 
flux in Lya is absorption by associated Hi. In §3.2.2, we 
saw that this Hi absorption will reduce the total flux of 
the Lya line by no more than 60%, which in most cases 
will not completely destroy the diagnostic value of the Lya 
flux for line ratios. In this section, we shall therefore first 
compare the observed ratios involving Lya with the values 
of the photo-ionization models found from the other UV 
lines, and then consider the geometrical and dust effects. 



Figure 14 presents two diagnostic diagrams with line- 
ratios of Lya, C IV, and C HI]. We find that the simple 
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Fig. 14. Line ratio diagnostic diagrams for Lya, C IV and C III]. Models are as in Figure |l^. Note the systematical 
over- luminosity of Lya at z > 3 (filled circles). 
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Fig. 15. Optical-UV line ratio diagnostic diagrams. Dots and lower limits are HzRG data, and open triangles are the 



Seyfert galaxies NGC 5643 and NGC 5728 (Evans et al. 1999). Models are as in Figure | 



photo-ionization model with a power law spectral index 
a = — 1 provides a good fit to the data in both diagrams. 
The range of the ionization parameter —2.5 <logio(U)< 
—2 is consistent with the values found from the dia- 
grams in Fig. |l^. Contrary to the diagrams in Fig. |l^, 
the a — —1.5 models do not provide a good fit to the 
data, because they overestimate the Lya flux. This might 
be due to geometrical effects: changing the viewing an- 
gle from front to back illuminated Lya can increase the 
C IV/Lya ratio by almost an order of magnitude (see Fig. 
7 of Villar-Martfn et al. (1996). We therefore do not take 



this better fit as evidence for a flatter spectral index of 
the incident ionizing continuum. 

Some individual HzRGs have highly discrepant 
Lya/C IV and Lya/C III] ratios. Both under- and over- 
luminous Lya occurs. Dust extinction has been proposed 
as an explanation for under-luminous Lya in the two most 
discrepant objects in our HzRG sample, TX 0211—122 
at 2; = 2. 
at z 



van Ojik et al. 1994 


) and MG J1019+0534 


2.765 ( 


Dey, Spinrad & Dickinson 1995, 



factor of ~3 lower compared to the best fitting photo- 
ionization models and the bulk of the other HzRGs. 
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However, in two objects with equally large amounts 
of dust, 4C 41.17 at ^ = 3.8 and TN J0121+1320 at 
z = 3.516, the Lya is a factor ~3 brighter than the models 
and the mean of the other HzRGs. TN J0205+2242 at 
z — 3.506 does not have a large dust mass detected, 
while it is also ^ 4x over- luminous in Lya. We conclude 
that there is no strong correlation between the detection 
of a large global dust mass and an anomalously low 
Lya flux. If dust obscuration is an important process in 
suppressing the Lya luminosity, it would have to be a 
localized process, but this would be hard to achieve, given 
the large extent of the Lya emission. High resolution 
Lya imaging, combined with sensitive imaging of the 
IR-dust emission (with ALMA) would be needed to 
examine this in more detail. Metallicity differences and 
geometrical effects probably play a more important role 
than dust obscuration in the destruction of Lya. The 
most remarkable observation from Figure |lj is that 
three of the seven HzRGs with over-luminous Lya in the 
C III]/Lya versus C IV/C HI] plot are at 2 > 3. We shall 
return to this in §6.3. 



5.5. UV-optical diagnostic diagrams 

In HzRGs, the rest-frame optical emission lines which 
are commonly used in line ratio diagnostic diagrams at 
low redshift, are redshifted into the near-IR. It has been 
a big technological challenge to obtain near-IR spectra 
of HzRGs with 2-4m telescopes, but recently, near-IR 
spectroscopy of HzRGs has become feasible with the ad- 
vent of efficient near-IR spectrographs on 8-lOm class 



ure is the exclusion of the pure shock models. Only 
4C 41.17 falls in a region of the plot that cannot be ex- 
plained by pure photo-ionization models, but the H/? de- 



telescopes (e.g., Larkin et al. 2000). It is to be expected 
that more near-IR spectroscopy of HzRGs will be avail- 
able soon. To date, a total of 15 z > 2 radio galaxies 
have published near-IR spectra obtained with the 3.8m 



UKIRT and 2.2m UH telescopes ( [E ales k RawHngs 1993 , 
Iwamuro et al. 1996| , Evans 1998 ), and recently the first 
near-IR sp ectrum obtained w ith NIRSPEC at Keck was 
pubhshed ( [Larkin et al. 20001 ). As shown by ADT98, a di- 
agnostic diagram of the optical [O III] to H/? versus the 
UV C IV/He II or C IV/C HI] line ratios can separate the 
shock and shock-|- precursor models. Such diagrams have 
the added advantage that they are relatively insensitive to 
the effects of dust extinction, as the ratios are determined 
from lines close in wavelength. In four objects (3C 256, 
MRC 2025-218, 4C 40.36 and 4C 41.17), H/3 was detected 
as well as [O HI], and two more objects (USS 0828-^193 
and 4C 48.48) have lower limits to the [O III] to H/3 ra- 
tio. The uncertainties in these points will be higher than in 
the optical spectroscopy because of the low signal-to-noise 
ratio and resolution of the near-IR spectroscopy. 

Figure |l^ reproduces the line-ratio diagrams from 
ADT98 with the five HzRG data points added. We 
also compare the HzRG data with combined HST and 
ground-based spec troscopy of two low-redshift Seyferts 
( Evans et al. 1999 ). The most firm conclusion from Fig- 



tection is marginal (Eales & Rawlings 1993). The other 
detection and the lower limits are consistent with the 
Seyfert points, and fall in a region where the simple 
photo-ionization models overlap with the highest values 
of the Ayi/i sequence and the high shock velocities of the 
shock -|-precursor models. A comparison with the range of 
parameter values of the UV diagnostic diagrams suggests 
that the simple photo-ionization models have more con- 
sistent values than the Am/i sequence, because the weak 
H/3, if confirmed by higher S/N observations, requires a 
higher fraction of matter-bounded clouds than the UV di- 
agrams suggest. Alternatively, the highest shock velocity 
(500 km s^^) predictions of the shock -f precursor models 
can also explain most of the points in Figure [l5|. This again 
would require an extension of these models to higher shock 
velocities to fully examine the importance of this ioniza- 
tion mechanism in HzRGs. 

To summarize, the comparison of the UV and optical 
line-ratios excludes the pure shock models, and is consis- 
tent with the high shock velocity shock -|-precursor models 
and the high values of the Am/i sequence. Pure photo- 
ionization models with a power-law spectral index a = — 1 
and ionization parameter U ~ 0.01 seem to provide the 
most consistent fit to HzRG spectra. Future near-IR spec- 
troscopy of HzRGs will be required to differentiate more 
unambiguously between these models. 

5.6. Nitrogen overabundance 

In the majority of HzRGs, the N V line was not de- 
tected. There are, however, several exceptions where N V 
is str ong, and even a few c ases where it is as bright as Lya 
(e.g., van Ojik et al. 1994 ). Neither shock ionization nor 
variations of the densities or ionization parameters in the 
photo-ionization models are able to produce such strong 
N V emission ( |Villar-Martm et al. 1999cD . The most likely 
explanation is that N V is over-abundant compared to the 
other species in HzRGs ( van Ojik et al. 1994 ). Fosbury et 
al. (1998, 1999) found that HzRGs follow a close correla- 
tion in a N V/He II vs. N V/C IV diagram which is par- 
allel to the relation defined by the broad line regions of 
quasars ( Hamann fc Ferland 1993| ) . These authors showed 
that this sequence can be explained by a large varia- 
tion of the metallicity (from Z ~ Zq to Z > IOZq) 
caused by rapid chemical evolution in a stellar popula- 
tion composed of massive stars. The parallel sequence de- 
fined by the extended emission line regions in HzRGs sug- 
gests that we also observe such a metallicity variation on 
much larger scales than the small central broad line region 
( [Fosbury et al. 1999] ). Vernet et al. (1999) found that the 



radio galaxy sequence is best modeled with a metallic- 
ity sequence in which the nitrogen abundance increases 
quadratically, while the other elements increase linearly. 
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Fig. 16. N V/ flux plotted C IV flux. The diagonal lines 
indicate constant N V/C IV ratios. Note that there is no 
strong clustering of the upper limits at high or low C IV 
fluxes, which would lead to an artificial correlation in the 
N V/C IV ratios. 



Figure presents the N V/He II vs. N V/C IV dia- 
gram with all the points from our HzRG sample, of which 
two thirds have only upper limits for N V. The inclusion 
of such a large number of upper limits could lead to red- 
shift selection effects. We believe this is a major problem 
in our sample, as N V is only observable over a limited 
redshift interval (z > 2), and the highest redshift objects 
have been observed longer or with larger aperture tele- 
scopes. This is also obvious in the flux-flux diagram of the 
N V and C IV lines, which shows a constant spread of the 
detections and upper limits on C IV flux. 

Using survival analysis, we find that the probability 
that the N V/C IV and N V/He II ratios are correlated 
with redshift is <35%, while the ratios are strongly corre- 
lated with each other (the generalized Spearman rank cor- 
relation coefficient is r = 0.95, signficance level > 99.99%). 
The extension of the radio galaxy sequence to both sides 
in Fig. |l^ suggests that the nitrogen abundance shows an 
even larger variation than the Z = OAZq to Z = 3Zq se- 
quence proposed by Vernct. The most extreme source in 
our sample is TXS J2353-0002 (De Breuck et al. 2000b), 
which has N V emission that is more than twice as lumi- 
nous as C IV and undetected He II emission, at least four 
times weaker than N V. 

Eventhough there is no statistical significance for a cor- 
relation between the N V/C IV ratio and redshift within 
the 2 < z < 4.4 redshift interval we can observe both these 
lines, it is clear from Figure p7| that none of the z > 3 radio 
galaxies occupy the high metallicity end of the diagram. 
Moreover, the limits from the VLT and Keck spectra on 



Fig. 17. N V/He II vs. N V/C IV. Radio galaxies at z < 3 
and z > 3 are represented with open and filled symbols, 
respectively. The dotted line represents the metallicit y 
sequence defined by quasars ( Hamann fc Fcrland 1993 ), 



with the numbers along the line representing the metal- 
licity in solar units. The solid line represents a metallicity 
sequence with quadratic N enhancement (U=0.035, 
power law spectral index a = —1.0) from Vernet et 
al. (1999). The sta rburst galaxy IRAS F10214-f4724 
( Serjeant et al. 199S ) is indicated with a cross and 
the hypcr-luminous, gravitationally lensed object 



SMM J02399-0136 ( |ViUar-Martfn et al. 1999cD is in- 
dicated by a star. Note that the z > 3 objects seem to 
fall preferentially in the low metallicity region of the 
sequence. 

the N V of TN J1338-1942 at z = 4.11 and 6C 0140-^326 
at z = 4.41 are quite strong (N V/C IV < 0.07), and 
suggest the z > 4 galaxies have even lower metalicities. It 
appears that the metalicities in the extended emission line 
regions of HzRGs show large variations, from values well 
below solar to several times solar. The highest metalicities 
only occur in z < 3 radio galaxies, while the most distant 
objects generally have sub-solar metalicities. 

6. Discussion 

6.1. Simultaneous shock and photo-ionization 

From the diagnostic diagrams of the UV line ratios, and 
the combined UV and optical line ratios in §5, we found 
that the line-ratios involving C IV, He II, and C III] can 
best be explained by photo-ionization models, whilst the 
C II] to C HI] ratio is better fitted by shock models with 
high shock velocities. Models that invoke a variation in the 
ratio of matter-bounded and ionization-bounded clouds 
provide a more consistent fit to all the UV line-ratios. 
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although the large scatter of the data around the model 
predictions does not appear to follow the Af^/i sequence. 

We find that within the ranges of emission line fluxes 
measured in our HzRG sample, neither redshift, nor ra- 
dio luminosity are important causes of this scatter. The 
radio size D also does not appear to influence the emis- 
sion lines, with the notable exception of the C II] /C III] 
ratio. This was interpreted by BRLOO in a C II] /C III] 
versus [Ne III] A 3869/ [Ne V] A 3426 diagram as evidence 
that the ionization mechanism in smaller radio sources is 
shocks, while the larger sources are photo-ionized. How- 
ever, the only object from their sample that also has a 
pubhshed C IV flux, 3C 324, lies in the shock-|-precursor 
region in their diagram, but cannot be explained by any 
of the shock models in our diagrams involving the C IV, 
He II, and C III] lines (Fig. H). The fluxes of these three 
lines in 3C 324 are very similar, putting the object con- 
sistently on the photo- ionization sequence with a — —1.0, 
and logj^Q U « —2.2 in all three plots involving these lines. 

We interpret this apparent inconsistency as an indica- 
tion that the integrated spectra of this, and most likely 
also other HzRGs, are a composite of differently ionized 
areas. The lower ionization state C II] line is much more 
sensitive to shock ionization than the other UV lines: in 
the shock ionization models of DS96, the C II] flux is sim- 
ilar to the C HI] flux, while in the pure photo-ionization 
models, C II] is ~ 5x weaker than C III]. An increasingly 
more dominant shock-ionized area will thus manifest itself 
first through brighter low-excitation lines such as C II]. 
We can now also understand why the A-^^i sequence pro- 
vides a much better fit to the Carbon-only plot than the 
pure photo-ionization models: from Table 2 of B96, we find 
that 25 X more C II] is produced in the matter-bounded 
clouds than in the ionization-bounded clouds; the addition 
of the ionization-bounded clouds could mimic the addition 
of a more dominant shock component. 

To examine this idea of a mixture of shock and photo- 
ionized gas in the integrated spectra in more detail, we 
have over-plotted a sequence of the fraction of shock 
to photo-ionization on the line-ratio diagrams in Fig- 
ure |l^. We picked the pure photo-ionization model with 
n — 100 cm~^, a ~ —1.0 and \ogU — —2.25, and added 
a contribution of the shock model with v = 400 km s^^ 
and Bj ^fn = 4yuG cm*^/^ in steps of 20%. This approach is 
likely to be a serious oversimplification, because we expect 
a range of parameters for both photo-ionization and shock 
models, but it might be a reasonable approach in a single 
object, such as 3C 324. A 40% shock contribution, in this 
object would provide an excellent fit in the Carbon-only 
plot, while it will not shift the predicted ratios in the other 
diagrams by as much. The C II] /C III] and C HI] /He II 
ratios show the largest variation when the fraction of shock 
ionization is increased, with the C II] flux increasing faster 
towards the shock models than the He II flux. A varying 
shock contribution might well explain some of the scat- 
ter observed in the diagrams of Figure [13^. Assuming that 



the a = —1.0 pure photo-ionization sequence is the most 
representative, the scatter in all diagrams always seems 
to fall along the side of the sequence pointing towards the 
(high shock velocity) shock models. Because the C II] line 
is the most sensitive to shock ionization, most points in 
the Carbon-only plot do not fall on the photo-ionization 
sequence, but closer to the shock models. 

Another effect which could raise the C II] /C HI] ra- 
tio is differential extinction. For a single, central photo- 
ionizing source (the AGN), the ionization parameter U of 
a cloud located closer to the source will be higher than for 
the same cloud located further out. If the higher excitation 
photo-ionized gas is more centrally located than the lower 
excitation gas, it will be more subject to extinction. How- 
ever, the C II] /C HI] ratio is not very sensitive to changes 
in the ionization parameter U for an ionizing continuum 
with an a = —1.0 power law spectrum. For example, to 
change the ionization parameter from logJ7 = —2.0 to 
logJ7 = —2.5 we need to move three times further out 
from the AGN, while the C II] /C HI] ratio only changes 
by a factor of two. We consider it very unlikely that this 
effect can be dominant, and so we neglect it. 

The above findings could undermine the diagnostic 
value of the different line-ratio diagrams. If the inte- 
grated HzRG spectra are a composite of different regions, 
we should not necessarily try to fit all lines with a sin- 
gle model. Photo-ionization remains the main process in 
HzRGs, but using only high excitation lines, we cannot 
exclude shock ionization, like VM97 and ADT98 did. The 
high excitation diagrams, which include lines that are in- 
sensitive to shock ionization, can be used to determine 
the best fitting parameters of the photo-ionization model, 
which can then subsequentially be used in the diagram 
which include lines that are more sensitive to shock ion- 
ization, such as C II] , to estimate the relative contribution 
of shock ionization. 

The study of BRLOO has shown that significant shock 
contributions only occur in sources with radio sizes 
< 150 kpc. Because only 15% of the z > 2 galaxies in 
our HzRG sample have radio sizes >150 kpc (and only 
6% >200 kpc), the contribution of shock ionization in 
HzRGs might well be important. There is also further 
kinematical and morphological evidence for the presence 
of shocks in HzRGs. At relatively low redshifts {z < 1), 
the line widths in radio galaxies are < 500 km s~^ 
and can also be explained by gravitational origins 
( [Baum fc McCarthy 2000| ). However, at z > 2, the line 
widths of the UV lines are ^ 1500 km s^^, which is more 
easily explained by acceleration due to the passage of a 
bow shock associated with the expansion of the radio 
source, although different processes might also play a 
role, such as (i) infall of material from large distances, 
(ii) large scale outflows of companion Lyman break 
galaxie s, or (iii) bipolar outflows produced b y super- 
winds ( Villar-Martfn, Binctte fc Fosbury 1999a). In some 
objects at low redshift (e.g., van Breugel et al. 1985 , 
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Fig. 18. Composite spectra from 3C sources at 



( |Bcst, Rottgcring fc Longair 2000a|), th e 3C and MRC 
samples (McCarthy & Lawrence 200C) and the MG 



(3tern et al. 1999). All spectra have been normalized to 
the C II] A 2326 line. The 3C galaxies in the top panel 
have the largest radio sizes, while the MG galaxies in the 
bottom panel the have smallest. This is clearly reflected 
in the dominance of the C II] line compared to the C IV, 
He II and C III] lines. 



Villar-Martm et al. 1999bD 
( iBickncll et al. 2000[ ), "there 
idence for the interaction of the 



and at high redshift 
is also morphological ev- 
radio jet with the 



ambient gas. 

The presence of a significant percentage of shock ion- 
ized gas in the integrated spectra of HzRGs would ex- 
plain the poor fit of the pure photo-ionization models to 
the C II] /C HI] ratio (Fig. |l3|). Because of the lack of 
sources with radio sizes >150 kpc, we do not find a pop- 
ulation of pure photo-ionization sources like BRLOO do. 
The radio size - ionization mechanism relation can also 
explain the different line ratios in three published com- 
posites of HzRGs (Fig. H). The MG sample of S99, which 
has an angular size cutoff d < 10", contains more sources 
with relatively strong C II] lines. The C II]/C HI] ratio 
in the MG composite spectrum of S99 is 1.6x stronger 
than the ratio in the 3C/MRC composite of McCarthy 
& Lawrence (2000), which does not have an explicit an- 
gular size cutoff. This can by easily interpreted if the 
MG sample contains more sources with important con- 
tributions of shock ionization. On the other end, a com- 
posite spectrum of 14 3CR galaxi es at 0.7 < z < 1.25 
( Best, Rottgering fc Longair 2000a ) contains more larger 
radio sources, and should therefore be the most represen- 
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Fig. 19. The Mg II/C HI] ratio plotted versus radio size. 
Note the absence of large sources with strong Mg II. 



tative composite of the photo-ionization in HzRGs, at least 
at z ~ 1. The C II]/C HI] ratio in this composite is 0.3, 
which is indeed much closer to the pure photo-ionization 
models. 

We can also use Fig. |l^ to identify other emission lines 
that could be sensitive to shock ionization. The only ob- 
vious candidate line is Mg II, which has a constant ratio 
compared to C II] in the 3C(z ^ 1) and 3C/MRC com- 
posites, and is very bright in the MG composite. How- 
ever, the high luminosity in the MG composite is most 
likely due to the high selection frequency of the MG sur- 
vey (5 GHz). At such frequencies, the sources will be more 
dominated by the flat spectrum radio core, and they will 
appear more like quasars, which have more prominent 
Mg II lines than radio galaxies, while t heir C II] lin es is 
five times weaker than Mg II (see e.g. 



Boyle 19901) ■ In- 



deed, S99 showed that the Mg II/C HI] ratio in the MG 
composite is similar to those in quasar composites, and 
deep spectro-polarimetric observations of radio galaxies 
and ultra-luminous IR galaxies have revealed a scattered 
broad Mg II component (Tran et al. 1998,2000). 

The plot of the Mg II/C III] ratio versus radio size D 
in our sample (Fig. |l^) looks remarkably similar to the 
C II] /C HI] versus D plot (Fig. The highest shock 
velocity models of DS96 indeed predict Mg II to be twice 
as strong as C HI] , while in a pure photo-ionization model 
with n = 100 cm"3, a = -1.0 and logC/ = -2.25, Mg II 
would be twice as faint as C HI]. We can therefore also 
use the Mg II line to identify a strong shock contribution 
to the integrated HzRG spectrum, but should be aware 
that this line is much more subject to small contributions 
from obscured broad line regions than C II] . 
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6.2. Emission-line luminosity - radio power correlations 6.3. The independent behavior of Lya 



In the previous section, we argued that the ionization 
mechanism in HzRGs is a composite of nuclear photo- 
ionization and shock ionization dominating the lower ion- 
ization lines. Both these processes are likely to be linked 
to the power of the central engine. To examine this de- 
pendence, we shall now consider the relation between the 
radio power and the emission line luminosity of the shock 
and photo-ionization sensitive lines. 

In §4.6.2, we found that the UV line luminosities are 
weakly correlated with the radio power over 2 orders 
of magnitude, with the strongest results for the C II] and 
Lya luminosity. This should be compared with the strong 
correlation between the Ha-f [Nil] or [O II] luminosity 
and radio power, which extends over nearly five orders of 
magnitude in radio power (e.g., Baum fc Hcckman 1989 



Rawlings fc Saunders 199lt [McCarthy 1993 , 

Willott et al. 1999|) . This relation has prompted these 



authors to suggest a common energy source (the AGN) for 
the total emission line luminosity and radio power. To cal- 
culate the total line luminosity, they have taken the flux 
in one or two of the brightest lines (Lya, [O II], [O III] or 
Ha) and calculated the flux in the other lines using fixed 
line ratios. In HzRGs, the brightest line is Lya, which can 
represent up to 1/3 of the total narrow line luminosity 
(assuming the line ratio from the composite of McCarthy 
1993). We can thus use Lya to calculate the total line 
luminosity, but this line has the disadvantage that is is 
highly influenced by resonant scattering and geometrical 
effects (e.g., Villar-Martin, Binette fc Fosbury 1996| ), 



which will cause the total line flux to be underestimated 
by a varying amount. The weaker correlation with 
the fainter high ionization lines can be explained by 
measurement errors and an intrinsic scatter caused by 
differences in the relation due to different accretion rates, 
environmental effects, or time variability of the AGN 
which reaches the narrow-line regions and radio lobes 
at different times (see Willot et al. 1999 for a detailed 
discussion) . 

The fact that the radio power does not correlate with 
the line ratios (Table ||), but does correlate with the 
individual UV lines, suggests that in the majority of 
HzRGs, the total emission line luminosity is dominated 
by one of the ionization processes (shocks or AGN photo- 
ionization), and this process has a common energy source 
with the radio power. Alternatively, it is possible that both 
mechanisms depend equally strong on radio power. 

It is remarkable that the shock ionization dominated 
C II] line luminosity has the strongest correlation with 
radio power (see Table |3|) . This suggests that in more lu- 
minous radio sources, the shock ionized regions become 
more prominent. However, the range in radio power in 
our HzRG sample is still limited, and deeper spectra of 
less radio-luminous radio galaxies are needed to confirm 
this trend. 



In §5.4, we found that radio galaxies with relatively over- 
luminous Lya emission occur more often at z > 3 than 
at lower redshifts. The correlation over the entire redshift 
range 2.0 < z < 4.4 is not very strong, and from Table |[ 
we find that it is not stronger than the correlation with 
radio size or radio power. However, the rest-frame Lya 
equivalent width also increases at z > 3. Figure |o| plots 
the relative intensity of Lya compared with the continuum 
(= equivalent width; this measure takes some account of 
absorption by dust), and with C IV. We find that at z > 
3, Lya is roughly twice as strong compared to both the 
continuum and C IV than at z < 3. We see similar trends 
in the ratios with He 11 and C 111], but we have too few 
data points to make a significant claim. This effect cannot 
be entirely due to selection effects, as HzRGs with weaker 
Lya could still have been detected (see Figure ^ . 

We interpret the over-luminous Lya as spectroscopic 
evidence for the youth of radio galaxies at z > 3. At these 
redshifts, the radio sources fo und from flux density-limite d 
surveys are inevitably young ( Blundell fc Rawlings 1999 ). 
It is therefore conceivable that they are still surrounded by 
large halos of primordial hydrogen from which the galaxy 
is formed and which trace the over-dense region that har- 
bours the AGN . The Lya halos at z > 3 are often very 
extended (e.g., van Ojik et al. 1996), so we may be wit- 



nessing the first epoch of star formation in primordial 
material. The very frequent occurrence of large amounts 
of H I absorption in the Lya profiles of the highest red- 
shift radio galaxies (§4.6.3) also indicates a larger amount 
of Hi surrounding the radio galaxy. As such, the higher 
hydrogen abundance at z ^ 3 indicates metallicity evolu- 
tion with redshift in forming massive galaxies. This would 
be consistent with our findings on the metallicity derived 
from the N V/C IV ratio (§5.6) and the recent results of 
Binette et al. (2000), who found strong evidence from co- 
spatial Lya and C IV associated absorption that the mate- 
rial in the outer halo must have a much lower metallicity 
{Z ~ O.OIZ0) than the more centrally located emitting 
material. 

At z < 3, the radio galaxies become more relaxed, and 
even the outer parts have experienced chemical enrich- 
ment through massive starbursts caused by the infall of 
the surrounding material on the central (proto-)galaxy, or 
by the passage of the radio jets. This picture is consistent 
with the change in rest-frame optical morphology seen in 
a sample of HzRGs observed in i^T— band by van Breugel 
et al. (1998). At z > 3, their morphologies are clumpy 
and radio-aligned, while at z < 3, they often have well 
organized elliptical profiles, much like their low redshift 
counterparts. 

7. Conclusions 

We have compiled a sample of 165 HzRG, half of which 
are at z > 2. From this sample, we have shown that: 
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Fig. 20. Lya equivalent width (left panel) and Lya/C IV ratio (right panel) as a function of rcdshift. Overplotted 
in each plot are three bins in redshift: z < 2.5, 2.5 < z < 3.5, and z > 3.5. The tickmarks on each bin are the 25 
percentile level, median and 75 percentile levels. Note that the z > 3.5 bin is always higher than the lower redshift 
bins, suggesting a higher Lya flux at the highest redshifts. 



• The different UV line luminosities and equivalent 
widths are strongly correlated, indicating that the nar- 
row line region gets most of its energy from a common 
energy source. The UV line luminosities also appear to 
be correlated with radio power, suggesting a common en- 
ergy source and either a dominating ionization mechanism 
or two mechanisms that increase equivalently with radio 
power. However, the shock ionization sensitive C II] line 
appears to be more strongly correlated than the photo- 
ionization dominated lines, suggesting the more powerful 
radio sources can ionize a larger amount of gas by means 
of shock ionization. 

• Using the systematic redshift determined from the 
He II line, we introduce a parameter Aj^y^ measuring the 
asymmetry in the Lya emission line due to Hi absorption. 
We find that the absorption occurs most frequently at the 
blue side, and mainly in sources with smaller radio sizes 
and at higher redshift. The higher Hi absorption could 
trace a denser surrounding medium which can confine the 
radio source more ( van Ojik et al. 1997 ), or it could reveal 
an un-pressurized, low density region that has not yet bee n 
affected by the passage of a radio jet ( Binettc et al. 2000| ). 

• We confirm the results from Fosbury et 
al. (1998,1999), who found that HzRGs occupy a 
sequence in a N V/C IV vs. N V/He II diagram which is 
parallel to the metallicity sequence defined by the BLR 
in quasars. The inclusion of upper limits to the N V flux 
further extends this sequence to much lower metallicities. 
The z > 3 objects occupy the lower metallicity part of 
this sequence. 

• For z > 3, Lya becomes over-luminous with respect 
to both the continuum and the other emission lines. We 



interpret this as indicating an increasing hydrogen abun- 
dance with redshift and evidence that the highest redshift 
radio galaxies tend to be in an earlier stage of their forma- 
tion epoch, still surrounded by the reservoir of primordial 
hydrogen from which the galaxy is forming. 

• Pure photo-ionization models with an ionizing con- 
tinuum of power law spectral index —1.0 provide better 
fits to the line ratios involving C IV, He II, and C III] than 
shock models, but cannot explain the C II]/C III] ratio. A 
combination of shock an photo-ionization in the integrated 
HzRG spectra can better explain the observed ratios, be- 
cause the C II] line is ~ 5 x more sensitive to shocks than 
photo-ionization. Such a shock ionization component will 
not be seen in high excitation line ratios, where the photo- 
ionization will dominate. Diagnostic line-ratio diagrams 
consisting of only high excitation lines will therefore fail 
to detect contributions of shock ionization. 

• The UV line ratios are not correlated with red- 
shift, radio size, or radio power, with one exception: 
high C II] /C III] ratios only occur in small radio sources 
(< 150 kpc). This confirms the results of BRLOO, who 
found that the high C II] /C III] ratios can be explained 
by shock ionization, but occur only in small radio sources, 
where the radio size is comparable to the extent of the 
emission lines regions. The C II] /C HI] ratio in HzRGs 
requires the addition of a component of shock ionized gas 
to the integrated HzRG spectra. This is consistent with 
the generally small radio sizes and high line widths ob- 
served. A composite spectrum of HzRGs with small radio 
sizes (MG) indeed shows a higher contribution of lines sen- 
sitive to shock ionization (C II], Mg II) than composites 
based on larger radio sources (3C, MRC). 



24 



Carlos De Breuck et al.: Statistical study of emission lines from HzRGs 



• Within the shock ionization models, the highest 
shock velocity models provide the best fit to the HzRG 
spectra. An expansion of these models to higher shock 
velocities and lower metallicities is needed to determine 
the full contribution of shock ionization to the integrated 
HzRG spectra. 
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Table A.l. High redshift radio galaxy sample 

source ~z LAS" Q'^ -P325) Log(Pi4oo) ^Lya*^ Line Flux Ref. 

" erg/s/Hz erg/s/Hz 10"^^erg/cm^/s 

















Lya 


N V 


C IV 


He II 


C III] 


C II] 




USS 0003-019 


1.541 


3.3 




35.23 


34.48 








5.90 


3.90 


3.40 


1.00 


DB 


BRL 0016-129 


1.589 


3.5 




36.09 


35.51 








1.60 




2.60 


1.20 


BRL 


MG 0018+0940 


1.586 


<1.2 




35.36 


34.84 








0.81 


0.42 


0.87 


0.65 


Ste 


MRC 0030-219 


2.168 


<0.3 




35.48 


34.93 




5.30 












Mcl 


6C 0032+412 


3.670 


3.0 




35.90 


35.16 




5.00 












Lac 


WN .10040+3857 


2.606 


1.4 




35.00 


34.14 


-0.13 


4.10 


<0.20 










DB 


MG 0046+1102 


1.813 


<1.2 




35.42 


34.79 








0.65 


0.55 


0.79 


0.74 


Stc 


3C 22.0 


0.935 


24.6 


1.2 


35.62 


35.04 












6.06 


1.25 


BOO 


BRL 0056-172 


1.019 


17.0 




35.70 


35.05 












1.40 


0.70 


BRL 


BRL 0101-128 


0.387 


16.0 




34.64 


34.17 
















BRL 


TN J0121 + 1320 


3.516 


0.3 




35.83 


34.97 
















DB 


MG 0122+1923 


1.595 


<1.2 




35.24 


34.77 








0.32 


0.38 


0.32 


0.23 


Stc 


BRL 0125-143 


0.372 


15.0 


1.7 


34.77 


34.24 
















BRL 


BRL 0128-264 


2.348 


33.0 




36.48 


35.81 




7.90 


<0.15 


1.10 




0.60 




BRL 


BRL 0132+079 


0.499 


9.0 




34.89 


34.43 
















BRL 


TXS ,70137+2521 


2.897 
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35.88 


35.16 
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MRC 0140-257 
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35.77 


35.13 




3.00 












Mc4 


6C 0140+326 


4.413 
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1.1 


36.03 


35.30 




6.40 


<0.01 


0.16 








Raw 


MRC 0152-209 


1.920 


1.0 




35.61 


35.03 




18.00 












Mc4 


MRC 0156-252 


2.016 


8.3 




35.61 


35.04 




15.00 
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USS 0200+015 


2.229 


5.1 




35.72 


34.97 
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35.01 
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6.70 
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USS 0211-122 
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1.8 


35.94 


35.21 


0.24 


5.70 


4.10 


5.60 


3.10 


2.20 




Rot 


USS 0214+183 
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1.6 
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Ste 
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35.79 
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20.00 
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3.70 


2.30 




Rot 
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USS 0417-181 
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3.7 


2.8 
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35.39 
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3.00 
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0.50 






Rot 


USS 0448+091 


2.037 


22.4 




35.31 
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12.20 
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1.20 


1.40 
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Rot 
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35.69 
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0.30 


0.10 


0.40 


DB 


4C 60.07 


3.788 


16.0 


2.8 


36.53 
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10.10 
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2.70 








Rot 
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0.357 


1.3 




33.33 
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DB 
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-0.05 


7.40 
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0.40 


0.60 


1.80 




Rot 
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3.153 
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35.41 


34.54 




0.80 


<0.10 










DB 


4C 41.17 
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20.0 


1.2 


36.45 


35.66 




14.60 


0.39 


1.32 


0.55 


0.91 




Rot 


B3 0731+438 


2.429 


10.8 


1.3 


36.16 


35.53 




31.00 












Mc3 


B3 0744+464 


2.926 


1.9 


1.4 


36.19 


35.56 




12.00 












Mc3 


USS 0748+134 


2.419 


13.2 




35.72 


34.97 




6.30 


<0.20 


1.80 


1.50 


1.40 




Rot 


WN J0813+4828 


1.274 


0.9 




34.42 


33.52 












0.20 


0.40 


DB 


USS 0828+193 


2.572 




1.2 


35.55 


34.84 


-0.06 


13.30 


<0.10 


1.90 


1.90 


2.00 




Rot 


BRL 0850-206 


1.337 


13.0 




36.00 


35.40 












2.10 


0.90 


BRL 


BRL 0851-142 


1.665 


7.0 




36.00 


35.47 








3.40 


2.30 


1.60 




BRL 


4C 12.32 


2.468 


15.0 










13.00 












GK 


USS 0857+036 


2.814 


4.0 




35.87 


35.16 


-0.44 


2.60 


<0.10 


1.00 


0.70 






Rot 


B2 0902+34 


3.395 


5.0 


1.2 


35.97 


35.41 




21.00 




2.00 








Lil 


3C 217.0 


0.898 


13.1 


3.3 


35.61 


34.99 












6.79 


5.34 


BOO 


TN J0920-0712 


2.760 


1.4 




36.00 


35.05 


-0.49 


44.00 


<0.10 


3.60 


3.00 






DB 


TN J0924-2201 


5.195 


1.2 




36.77 


35.74 




0.35 


<0.01 










DB 


TN J0941-1628 


1.644 


1.9 




35.71 


34.87 








3.20 


0.90 


2.00 


1.90 


DB 


3C 226.0 


0.818 


31.8 


1.3 


35.60 


34.96 












1.94 


1.00 


BOO 


USS 0943-242 


2.923 


3.7 




36.10 


35.36 


-0.37 


20.10 


1.70 


3.90 


2.70 


2.30 




Rot 


MG 1019+0534 


2.765 


2.2 


1.3 


35.66 


35.19 


-0.68 


0.84 


0.23 


1.04 


0.85 


0.49 




Ste 


TN ,71033-1339 


2.427 


2.0 




35.89 


35.02 


0.17 


9.80 


<0.10 


2.30 


0.80 


0.70 




DB 


BRL 1039+029 


0.535 


6.4 




34.98 


34.56 
















BRL 


3C 247.0 


0.749 


13.9 


1.5 


35.35 


34.89 












2.73 


1.18 


BOO 


TN J1102-1651 


2.111 


3.0 




35.57 


34.72 


-0.78 


2.70 


<0.20 


1.00 


1.30 


1.10 




DB 


3C 252.0 


1.104 


56.7 


2.0 


35.74 


35.04 










6.16 


3.56 


1.38 


BOO 


TN ,71112-2948 


3.090 


9.1 


1.5 


36.00 


35.12 


0.22 


2.90 


<0.05 


0.70 


1.20 






DB 


USS 1113-178 


2.239 


10.3 


2.1 


35.72 


35.08 




6.40 


<0.20 


1.70 


0.70 


2.80 




Rot 


3C 256.0 


1.824 


4.0 


2.9 


36.15 


35.51 




54.20 


1.40 


5.23 


5.47 


4.28 


2.03 


Sim 


WN J1123+3141 


3.217 


25.8 


1.6 


35.99 


35.06 


-0.18 


6.20 


<1.00 


3.40 


2.40 






DB 


USS 1138-262 


2.156 


11.4 


1.5 


36.34 


35.57 




13.90 


<0.30 


0.80 


1.30 


1.30 




Rot 


BRL 1138+015 


0.443 


5.2 




34.80 


34.39 
















BRL 


MG 1142+1338 


1.279 


<0.4 




34.88 


34.50 
















Ste 


BRL 1140-114 


1.935 


3.9 




36.30 


35.60 








1.00 


0.50 


0.60 


1.90 


BRL 


B2 1141 + 354 


1.781 


11.0 




35.45 


34.82 




47.00 








0.80 




AS 


3C 2()5.0 


0.810 


78.8 


1.5 


35.80 


35.07 












22.78 


3.56 


BOO 


4C 26.38 


2.608 


23.0 














8.90 


5.70 


2.40 




Rot 


B2 1230+349 


1.533 


11.0 




35.44 


34.85 












3.00 


2.00 


AS 


6C 1232+39 


3.220 


8.0 




36.28 


35.48 




10.00 












Lac 


VLA ,7123642+621331 


4.424 


0.2 










0.07 












Wad 


USS 1243+036 


3.570 


7.0 


1.4 








23.50 


<0.50 










Rot 
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Source 




z 


LAS° 






Log(P325) 

org/s/Hz 


Log(Pi40o) 

crg/s/Hz 




Lya 


N V 


Line Flux 
10^^'^org/cm^/s 
C IV Ho II 


C III] 


C II] 


Rcf. 


WN J1333+3037 


1 


213 





4 






34.94 


33.65 

















30 


0.40 


DB 


WN J1338+3532 


2 


769 


11 


6 






35.72 


34.97 


-0.01 


17.80 


<0.60 


1 


30 


3.00 


2 


20 




DB 


TN J1338-1942 


4 


110 


1 


4 


2 


7 


36.29 


35.46 


-0.50 


10.10 


<0.07 


1 


00 


0.20 








DB 


3C 289.0 





967 


10 


6 


1 


1 


35.65 


35.09 














2 


81 


1.01 


BOO 


BRL 1344-078 





384 


<10 









34.72 


34.15 




















BRL 


4C 24.28 


2 


879 


7 





1 


1 


36.39 


35.65 




7.30 


6.90 


1 


70 










Rot 


WN J1356+3929 





253 


1 


6 






33.67 


32.86 




















DB 


BRL 1354+013 





819 


33 































BRL 


USS 1357+007 


2 


673 


4 


6 










-0.06 


9.10 


<0.20 


1 


70 










Rot 


MG 1401+0921 


2 


093 


3 


5 






35.53 


34.92 











41 


0.50 





34 


0.17 


Sto 


TN J1402-1510 





739 


15 


9 






35.00 


34.38 




















DB 


3C 294.0 


1 


786 


14 


5 


1 


4 








155 


3.10 


15 


50 


15.50 


18 


60 




Mc2 


USS 1410-001 


2 


363 


23 


5 


1 





35.65 


35.00 


0.06 


28.00 


1.68 


3 


36 


2.52 


1 


12 


0.36 


Cim 


BRL 1411-057 


1 


094 


47 





1 


4 


35.72 


35.06 














1 


10 




BRL 


BRL 1422-297 


1 


632 


<10 









36.11 


35.57 








4 


30 


2.10 


1 


00 




BRL 


USS 1425-148 


2 


349 


10 


7 






35.88 


35.28 


-0.01 


20.70 


<0.70 


2 


30 


2.30 


1 


00 




DB 


8C 1435+635 


4 


250 


4 


3 


1 


1 


36.95 


36.12 




1.50 
















Spi 


USS 1436+157 


2 


538 


4 


7 


1 


5 






-0.01 


42.00 


<0.80 


17 


00 


6.00 


9 


40 




Rot 


BRL 1436-167 





146 


<12 









33.77 


33.29 




















BRL 


BRL 1509+015 





792 


7 


2 


1 


4 


35.37 


34.87 




















BRL 


USS 1545-234 


2 


755 


6 


3 


1 


9 


35.92 


35.24 


0.01 


6.70 


<0.20 


1 


10 


1.80 








Rot 


3C 324.0 


1 


208 


11 


1 


1 


7 


35.92 


35.34 








3 


67 


2.70 


3 


47 


1.43 


BOO 


USS 1558-003 


2 


527 


7 


7 










0.18 


14.90 


<0.10 


2 


70 


1.70 


1 


20 




Rot 


BRL 1602-174 


2 


043 


37 





3 


2 


36.32 


35.67 








10 


00 


4.80 


2 


70 




BRL 


BRL 1602-288 





482 


61 





1 





35.09 


34.48 




















BRL 


BRL 1602-093 





109 


290 































BRL 


BRL 1603+001 





059 


11 









32.96 


32.22 




















BRL 


BRL 1621-115 





375 


<20 









34.73 


34.25 




















BRL 


3C 340.0 





775 


44 


7 


1 




35.47 


34.84 














7 


08 


0.75 


BOO 


BRL 1628-268 





166 


93 































BRL 


BRL 1643+022 





095 


7 


1 


1 


























BRL 


TXS J1650+0955 


2 


510 


19 


4 










0.28 


20.90 


<0.10 


3 


20 


2.70 


1 


20 




DB 


BRL 1649-062 





236 


85 































BRL 


WN J1703+3739 





256 


30 


1 






31.98 


31.41 




















DB 


USS 1707+105 


2 


349 


21 


7 












4.30 


<0.10 






0.90 





20 




Rot 


3C 352.0 





806 


12 


3 


1 




35.50 


34.84 














4 


11 


2.36 


BOO 


3C 356.0 


1 


079 


73 


4 






35.72 


35.04 












4.33 


4 


41 


0.52 


BOO 


BRL 1732-092 





317 


45 









34.44 


34.00 




















BRL 


8C 1736+650 


2 


400 


17 









35.08 


34.24 




18.00 
















L99 


7C 1740+6640 


2 


100 


<0 


5 






34.84 


34.50 




4.10 







30 










L99 


7C 1758+6719 


2 


700 


45 









35.40 


34.69 




12.00 













80 




L99 


7C 1802+6456 


2 


110 


26 









35.54 


34.76 




7.10 




1 


00 







80 




L99 


8C 1803+661 


1 


610 


36 









34.55 


33.93 








5 


30 


2.60 


1 


90 




L99 


3C 368.0 


1 


132 


8 


5 


2 




35.93 


35.10 












5.11 


10 


57 


9.98 


BOO 


7C 1805+6332 


1 


840 


14 









35.13 


34.47 




17.00 




2 


50 










L99 


7C 1807+6719 


2 


780 


1 


9 






35.31 


34.78 




2.10 


0.70 














L99 


4C 40.36 


2 


265 


20 









36.24 


35.42 








6 


20 


5.60 


5 


90 




Rot 


BRL 1859-235 


1 


430 


4 


2 






36.29 


35.68 








3 


40 


4.60 


4 


70 


6.70 


BRL 


TXS J1908+7220 


3 


530 


14 


4 






36.44 


35.64 




32.00 


<0.20 














DB 


WN J1911+6342 


3 


590 


1 


8 






35.56 


34.66 


-0.64 


1.40 


<0.03 





08 


0.08 








DB 


BRL 1912-269 





226 


48 





1 




34.33 


33.57 




















BRL 


BRL 1920-077 





648 


23 









35.19 


34.62 




















BRL 


4C 48.48 


2 


343 


17 





1 




35.87 


35.20 








6 


10 


3.70 


2 


80 




Rot 


MRC 2025-218 


2 


630 


5 


1 






35.89 


35.25 




4.00 


0.62 





69 


0.35 





97 




VM99 


TXS J2036+0256 


2 


130 


3 


1 






35.86 


35.09 


0.01 


6.80 


<0.30 





60 


0.70 


1 


20 




DB 


MG 2037-0011 


1 


512 


<0 


4 






34.93 


34.63 

















44 


0.10 


Sto 


MG 2058+0542 


1 


381 


<0 


4 






35.40 


35.02 

















71 


0.41 


Ste 


MRC 2104-242 


2 


491 


21 


8 






36.26 


35.41 




57.00 


<3.80 


3 


80 


1.90 


2 


66 




VM99 


4C 23.56 


2 


483 


47 





1 










8.00 


1.36 


2 


08 


1.52 


1 


28 


0.80 


Cim 


MG 2109+0326 


1 


634 


<1 


2 






35.28 


34.71 












0.32 





21 


0.23 


Ste 


MRC 2115-253 


1 


114 


1 


7 






35.16 


34.45 














37 


00 




Mcl 


MG 2121 + 1839 


1 


860 


6 


3 






35.32 


34.73 











53 


0.14 





24 


0.15 


Ste 


BRL 2120-166 





882 


14 









35.49 


34.86 


















1.60 


BRL 


MG 2144+1928 


3 


592 


8 


5 






36.37 


35.65 


-0.47 


6.15 


<0.04 





58 


0.35 








Ste 


USS 2202+128 


2 


706 


3 


4 


1 


4 


35.79 


35.10 




7.70 


<0.10 


1 


90 


2.40 








Rot 


30 441.0 





708 


26 





3 


1 


35.38 


34.82 


















0.30 


BOO 


MRC 2224-273 


1 


679 





4 






35.37 


34.72 




47.00 
















Mc4 


MRC 2247-232 


1 


326 


9 


3 






35.57 


35.01 














1 


50 




Mc4 


USS 2251-089 


1 


986 


4 


5 


1 




35.66 


35.07 








3 


30 


1.30 


1 


50 




Rot 


TN J2254+1857 


2 


153 


2 


7 






34.97 


34.15 




33.00 


<1.00 






3.20 








DB 


MRC 2303-253 





740 


18 


6 






35.05 


34.35 














7 


00 




Mcl 


MG 2308+0336 


2 


457 










35.89 


35.32 


-0.74 


2.93 


0.57 





63 


0.39 





45 


0.83 


Sto 


BRL 2318-166 


1 


414 


<5 









36.20 


35.52 












7.80 


6 


20 


2.00 


BRL 


MRC 2318-244 


1 


113 


25 





1 




35.29 


34.74 














7 


00 




Mc4 


TXS J2321+2237 


2 


553 


8 


1 






35.90 


35.14 


0.17 


8.10 


<0.20 














DB 


BRL 2322-052 


1 


188 


7 


8 


























1.80 


BRL 


TXS J2334+1545 


2 


480 


6 


4 






35.69 


34.98 




3.10 


<0.10 





10 







80 




DB 


BRL 2347-026 


1 


036 


<2 









35.56 


35.03 


















0.60 


BRL 


TXS J2351 + 1034 


1 


334 


2 


2 






35.39 


34.60 




















DB 


4C 28.58 


2 


891 


28 





2 


3 


36.30 


35.44 











30 


1.60 


1 


80 




Rot 


TXS J2355-0002 


2 


587 


33 


8 


1 


6 


35.96 


35.28 




1.30 


0.80 





30 










DB 



°' Radio largest angular size, see §3.1.2 
^ Radio lobe distance ratio, see §3.1.3 

Lya asymmetry, see §3.2-2 

REFERENCESi AS-|lAllington Smith ot al. 
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Table B.l: Spearman rank correlation coefficients and significance levels 

' Q ^325 -^Lyo -^Lyo -^CIV 'n.U ^CIII] "'"yt ""ciV "^HeH ""ciQI ^"CIV ^"Hell ^"CIII] 





















> 






























146 


***** 


47.60 


5 9.99 


99.97 


99.97 


79.24 


88.00 


98.95 




34.31 


52 


,49 


57.99 


31 .69 




89.02 


46 


,1 1 


39 


,83 


33. 15 


1 3 


,72 


0.09 
54 


-0.09 

53 





94.83 


87,89 


34.85 


65.33 


8. 15 


74,68 




68.68 


96 


,79 


22.26 


99. 10 




77.43 


93 


,96 


6 


,82 


14.59 


72 


,01 


0.65 
147 


-0.08 
128 


0.28 
49 





> 99.99 


9.88 


99.63 


99. SS 


99,47 


> 


99.99 


70 


,25 


69.60 


91.60 




93.39 


42 


,68 


50 


,97 


0.23 


30 


,43 


-0.55 
147 


0.32 
128 


-0.22 
49 


-0.41 


***** 


0.87 


59.86 


97.70 


98,39 




99.93 


41 


,11 


99.24 


99,78 




99.47 


77 


,35 


56 


,62 


58.46 


30 


,43 


-0.34 

30 


0.65 
27 


-0.16 
10 


-0.03 

25 


-0,00 

25 


***** 


70.02 


83.46 


91,37 




18.12 


56 


,33 


2.25 


94.43 




73.50 


89 


,20 


82 


,18 


71.79 


72 


,51 


0.02 
80 


0.15 
74 


-0.19 

27 


0.35 
67 


-0, 10 
67 


0.20 
30 


***** 


> 99.99 


> 99,99 


> 


99.99 


99. 


19 


19. 13 


42.28 




79.98 


96 


,28 


91 


,73 


76.35 


55 


,71 


0.26 
73 


0.20 
63 


0.02 
28 


0.40 
64 


— 0.28 
64 


0.27 
27 


0.77 
50 


***** 


> 99,99 


> 


99.99 


28 


,92 


38.24 


54.77 




81.86 


88 


,07 


70 


,85 


7.87 


29 


,53 


0.30 
70 


0.33 
59 


-0.23 
27 


0.35 
61 


— 0.31 
61 


0.35 
25 


0.65 
40 


0.83 
60 




> 


99.99 


36 


,96 


81.67 


81.67 




99.19 


82 


,35 


2 


,98 


73.86 


16 


,81 


0.43 
88 


0.05 
74 


0.19 
31 


0.47 
79 


-0,37 
79 


0.06 
17 


0.64 
35 


0.73 
55 


0.83 
59 






10 


,78 


97.49 


91.50 




94.78 


82 


,05 


44 


,46 


38.27 


83 


,29 


0.11 
49 


-0.11 

45 


-0.46 

22 


0. 17 
41 


-0,09 
41 


-0.18 

20 


0.38 
47 


0.07 
31 


0, 10 
24 




-0.03 
19 






99.96 


99,99 




94.38 


35 


,50 


10 


,82 


43.30 


84 


,15 


-0.21 
55 


-0.12 
47 


-0.07 
21 


-0. 15 
49 


0,38 
49 


0.01 
24 


0.04 

35 


-0.07 
53 


— 0,20 
46 




-0.35 
41 





,65 
25 




> 99,99 


> 


99.99 


89 


,10 


10 


,84 


83.94 


95 


,98 


-0.05 
60 


-0.06 
50 


-0.53 
23 


-0.24 
53 


0,41 

53 


-0.42 
21 


-0.11 
28 


-0. 11 
49 


— 0,09 
58 




-0.24 
51 





,73 
22 


0.78 
45 




> 


99.99 


99. 


45 


94 


,84 


67.66 


91 


,39 


-0.15 
75 


-0.21 

62 


-0.23 
30 


-0.22 
69 


0,33 
69 


-0.31 
1 5 


-0.27 
24 


-0.20 
45 


— 0,37 
50 




-0.23 
75 





,50 

1 5 


0.80 
38 


0.84 
48 






99. 


76 


92 


,21 


96.29 


76 


,07 


-0.26 

55 


0.09 
51 


0.43 
20 


-0.09 
45 


— 0,18 
45 


0.30 
29 


0.29 
53 


0.25 
41 


0.23 
35 




0.27 
27 


-0, 


08 

35 


-0.28 

33 


— 0.52 
27 




-0.64 

20 






99. 


73 


92.78 


88 


,84 


-0.26 
56 


-0.08 
5 


0.02 
1 9 


0. 1 

5 


0. 11 

50 


0.28 
25 


0.28 
40 


0. 15 
5 4 


— 0.01 
48 




0.1 
40 


-0, 


03 
27 


0.02 
44 


— 0.31 
40 




-0.31 

33 





,47 
38 






90.13 


99 


,15 


-0.03 
50 


-0.07 
44 


-0.05 
16 


-0.00 
43 


-0, 13 
43 


0.23 
23 


0.21 
34 


0.02 
45 


0, 17 
48 




0.08 
38 





,13 
21 


-0.24 
36 


-0. 16 
40 




-0.38 
31 





,32 
33 





,26 
43 




16 


,92 


0.01 
47 


0.03 
38 


-0.32 
13 


0.06 
40 


-0,20 
40 


-0.28 
17 


-0.15 

28 


0.06 
39 


0,04 
39 




0.21 
46 


-0, 


36 
17 


-0.38 

29 


-0.30 

33 




-0.20 
38 





,31 
27 





,44 

35 


0. 17 
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